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ABSTRACT 


Over  the  last  two  decades,  several  nonconventional  gyro  schemes  have 

been  proposed.  Directed  primarily  at  eliminating  the  bearing  problems 

in  conventional  mechanical  gyros,  it  was  anticipated  that  either  better 

performance  could  be  achieved  or  performance  comparable  to  conventional 

gyros  could  be  achieved  at  a  lower  cost.  No  attempt  is  made  here  to 

compare  or  contrast  all  of  these  various  schemes,  but  rather  to  look 

3 

at  certain  aspects  of  a  particular  device,  a  cryogenic  He  nuclear 
gyro,  and  to  provide  models  for  estimating  its  performance. 

The  nuclear  gyro  is  highly  dependent  upon  the  stability  and  homo¬ 
geneity  of  the  magnetic  field  environment.  A  combination  of  low  temper¬ 
ature  physics  technologies  makes  such  a  device  practical.  By  generating 
the  applied  field  with  superconducting  coils  the  stability  of  this 
field  is  guaranteed  by  the  inherent  stability  of  the  persistent  super¬ 
currents  in  these  coils.  Then  protection  from  a  changing  external  en¬ 
vironment  is  provided  by  enclosing  the  device  within  a  superconducting 
shield . 


Nevertheless,  even  in  an  ideal  magnetic  field  environment,  the 
nuclear  gyro  is  subject  to  various  rectification  mechanisms.  Of  par¬ 
ticular  interest  is  the  role  kinematic  rectification  plays  in  this 
device.  Rectification  mechanisms  are  discussed,  and  mathematical 
models  derived  for  calculating  the  methods  of  kinematic  rectification 


drift.  An  experiment  designed  to  measure  and  distinguish  the  kine¬ 
matic  rectification  effects  are  discussed.  /v  -X 
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Chapter  I 


INTRODUCTION 


A.  BACKGROUND 

Bearing  friction  has  long  been  recognized  as  a  limiting  factor  on 
lifetime  and  accuracy  for  conventional  mechanical  gyroscopes.  The  evo¬ 
lution  of  these  devices  has  proceeded  through  many  improvements  in 
ball  bearings,  raceways,  and  lubricants,  then  gas  bearing  technologies, 
and  more  recently,  electrostatic  and  electromagnetic  suspension  schemes. 
The  gas  bearing,  as  well  as  the  electrostatic  and  electromagnetic 
suspensions  employ  very  small  gaps  between  the  rotating  member  and  its 
support  housing.  Contaminants  in  these  gaps,  which  may  originate  during 
assembly,  from  outgassing  materials  in  the  vicinity  of  the  gap,  or  from 
fretting  out  of  particulate  matter  during  contact  from  starts,  stops, 
or  over-slewing,  degrade  performance  and  cause  life  problems. 

Two  nonconventional  schemes  were  proposed  in  the  early  1960's  to 
eliminate  bearing  problems  by  eliminating  the  need  to  support  a  spinning 
mass.  These  were  the  ring  laser  gyro  and  the  nuclear  magnetic  resonance 
gyro.  There  are  no  rotating  parts,  no  gimbals  and  thus  no  bearings  in 
either  of  these  types.  This  aspect  alone  probably  provided  sufficient 
impetus  for  development;  but  there  was  also  the  potential  for  faster 
warm-up  times,  lower  costs,  and  the  ability  to  survive  more  severe  shock 
and  vibration  environments  than  obtainable  from  mechanical  gyros  of 
comparable  accuracy. 

To  date  the  ring  laser  gyro  has  already  realized  much  of  this  poten¬ 
tial,  with  the  Honeywell  version  now  operational  on  the  Boeing  767  air¬ 
craft.  The  basic  concept  [Hacek  1963]  makes  use  of  the  high  sensitivity 
of  the  oscillating  frequency  of  a  laser  to  dimensional  variations  of 
its  resonant  structure.  The  particular  structure  employed  is  that  of 
a  closed  optical  path  (ring)  in  which  light  may  travel  in  either  direc¬ 
tion  (Fig.  1-1).  With  no  rotation  of  the  ring,  the  clockwise  and 


-1- 


counterclockwise  beams  travel  the  same  path  length  in  completing  a  pass 
through  the  ring.  However,  for  rotations  about  an  axis  normal  to  the 
plane  of  the  ring,  the  effective  path  lengths  per  pass  are  different 
for  the  two  beams.  This  change  in  the  effective  lengths  of  the  resonant 
structures  leads  to  different  resonant  frequencies  for  the  clockwise  and 
counterclockwise  beams.  Then  to  first  order  this  frequency  difference 
is  proportional  to  the  component  of  the  rotation  normal  to  the  plane  of 
the  ring. 


FIG.  1-1  BASIC  ELEMENTS  OF  A  RING  LASER  GYRO 


The  basic  concept  of  the  nuclear  gyro  is  more  closely  related  to 
the  mechanical  gyro.  Certain  nuclear  species  possess  intrinsic  spin 
angular  momentum;  and  this  angular  momentum  can  serve  as  an  inertial 
reference  in  the  same  manner  that  the  wheel  angular  momentum  does  for 
a  mechanical  gyro.  However,  since  the  nucleus  cannot  be  gimballed, 
gimbal  angles  are  not  used  to  read  out  the  spin  axis  direction.  Instead, 
advantage  is  taken  of  the  fact  that  these  nuclear  species  possess  a 
magnetic  dipole  moment  which  is  always  parallel  (or  always  antiparallel) 
to  the  spin  direction.  Thus  a  sensitive  magnetometer  can  be  used  to 
determine  the  spin  direction. 

3 

The  particular  device  under  investigation  here,  the  cryogenic  He 

nuclear  gyro,  had  yet  another  raison  d'etre,  in  addition  to  the  previously 

mentioned  potential  advantages.  Physicists  have  long  searched  for  an 

electric  dipole  moment  in  neutrons.  It  was  proposed  [Fairbank  1967]  that 

3 

a  nuclear  gyro,  using  He  as  the  nuclear  species,  with  the  addition  of 
an  electric  field,  would  provide  an  ideal  instrument  for  the  dipole 

3 

moment  search  (the  He  nucleus  contains  a  single  unpaired  neutron).  This 

3 

proposal  constituted  the  motivation  for  the  design  of  a  cryogenic  He 
gyro  to  accomplish  this  end  [Taber  1978];  and  Taber's  work  then  con¬ 
stituted  the  motivation  for  this  thesis. 

B_,  BACKGROUND  ON  PREVIOUS  RESEARCH 

The  Singer  Company,  Kearfott  Division,  began  work  on  a  nuclear 

gyro  around  1960.  Their  device  [Simpson  19793  uses  two  stable,  odd 

199  201 

isotopes  of  mercury  (  Hg  -  Hg)  as  nuclear  species  and  a  rubidium 

vapor  optical  pumping  magnetometer  to  sense  the  spin  axis  directions  of 

the  two  nuclear  species.  The  best  published  results  to  date  show  a  rate 

bias  stability  of  +_  0.312  deg/h.  The  Hg  species,  which  has  a  nuclear 

spin  of  3/2  fi  and  thus  possesses  a  nuclear  quadrupole  moment.  EJectric 

field  gradients,  present  at  the  walls  of  the  sample  cell  or  resulting 

from  the  light  beams  used  for  optical  pumping,  interact  with  this  nuclear 

201 

quadrupole  moment  causing  the  Hg  spin  axes  to  precess.  This  is  thought 
to  be  the  major  perturbation  to  the  rate  bias  stability  in  this  device. 
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Litton  Guidance  k  Control  Systems  began  work  on  a  similar  device 

129 

around  1968,  Their  device  [Kanegsberg  1979]  uses  two  noble  gases  (  Xe  - 
83 

Kr)  as  nuclear  species  and  also  employs  a  rubidium  vapor  optical 

199  129 

pumping  magnetometer.  Like  Hg,  Xe  has  a  nuclear  spin  of  l/2  & 

83 

and  thus  possesses  no  nuclear  quadrupole  moment.  Kr,  however,  is  a 
spin  9/2  Ti  species  thus  subject  to  nuclear  quadrupole  interactions 
with  electric  field  gradients.  However,  a  more  serious  problem  found 
with  this  device  is  due  to  a  hyperfine  interaction  of  the  rubidium 
atoms  with  the  noble  gas  atoms.  The  result  of  this  interaction  is  the 
presence  of  an  effective  magnetic  field  at  the  noble  gas  nucleus  during 
collision  with  a  rubidium  atom.  Since  the  magnitude  of  this  effective 
field  is  different  for  the  two  nuclear  species,  a  rate  bias  results. 

This  interaction  was  also  found  to  cause  large  sensitivities  to  optical 
pumping  lamp  intensities  and  to  temperature.  Kanegsberg' s  report  does 
not  present  any  rate  bias  stability  data  but  does  indicate  that  the 
present  model  can  be  easily  compensated  to  yield  10  deg/h, 

3 

The  cryogenic  He  nuclear  gyro,  being  developed  here  at  Stanford 

3  3 

University,  employs  only  a  single  species,  He  ,  Since  He  is  a  spin 
l/2  tl  species,  it  possesses  no  nuclear  quadrupole  moment.  Then,  since 
the  rubidium  vapor  optical  pumping  magnetometer  is  not  used,  this  device 
is  not  subject  to  the  hyperfine  interaction  mentioned  previously.  The 
major  difficulty  with  any  single  species  device  is  that  it  makes  no 
distinction  between  a  shift  in  the  ambient  magnetic  field  and  an  actual 
rotation.  Thus  it  is  imperative  to  strictly  control  the  magnetic  environ¬ 
ment.  Fortunately  much  successful  work  had  preceded  this  effort  to 

establish  very  low,  homogeneous,  stable  magnetic  fields  inside  super- 

3 

conducting  shields  [Cabrera  1975],  Another  difficulty  with  the  He 

nuclear  gyro  is  that,  at  the  cryogenic  operating  temperatures,  optical 

pumping  is  not  possible.  Thus  long  relaxation  times  (see  Ch.  II)  are 

3 

necessary  for  a  practical  He  nuclear  gyro.  A  longitudinal  relaxation 
time  of  140  hours  was  obtained  in  a  pyrex  sample  cell  under  the  cryogenic 
environment  proposed  for  this  gyro  [Taber  1978]  .  Taber  has  also  per¬ 
formed  the  analysis  to  predict  the  effects  that  dimensional  variations, 
temperature  variations,  magnetometer  noise,  gyro  housing  asymmetries,  and 
London  moment  will  have  on  the  bias  stability  of  this  gyro. 


I 


£.  STATEMENT  OF  THE  PROBLEM 

The  purpose  of  this  research  is  to  analyze  the  expected  performance 
3 

of  the  He  gyro  with  respect  to  vibrations,  rotations,  and  accelerations, 

such  as  would  be  encountered  in  typical  gyro  applications.  Previous  work 

3 

on  the  cryogenic  He  nuclear  gyro  involved  a  detailed  development  of  the 

theory  of  relaxation,  relaxation  measurements  under  a  number  of  situa- 

3 

tions,  and  an  analysis  of  the  expected  performance  of  the  He  gyro  with 

3 

respect  to  the  requirements  of  the  He  electric  dipole  moment  experiment 
[Taber  1978],  For  the  purposes  of  that  experiment,  only  a  single  gyro — 
with  its  input  axis  aligned  with  the  earth's  rotational  velocity — was 
needed.  Mechanical  motions  were  deliberately  kept  to  a  minimum,  i.e., 
to  the  extent  that  the  experiment  was  conducted  on  a  concrete  isolation 
pad  which  was  separate  from  the  building  floor. 

D.  SUMMARY  OF  PRESENT  WORK 

The  principal  contributions  of  this  research  are  listed  below, 

•  A  theoretical  analysis  of  the  effects  of  cross  axis  inputs,  in¬ 
cluding  a  prediction  that  certain  oscillatory  rate  inputs  could 
precess  the  sample  magnetization  toward  the  applied  magnetic 
field,  thus  reducing  the  signal  to  noise  of  the  gyro. 

•  A  theoretical  analysis  of  rectification  of  oscillatory  inputs 
including  both  rectification  of  oscillatory  cross  axis  inputs 
and  kinematic  rectification. 

•  Modification  of  Taber's  laboratory  equipment  provided  a  device 
which  could  be  tested  as  a  nuclear  gyro. 

•  Design  of  an  experiment  to  use  the  above  device  to  verify  the  pre¬ 
dicted  kinematic  rectification  effects  which  produce  drift  in 

the  presence  of  coning  motion. 

•  A  theoretical  analysis  of  the  effects  of  acceleration  and  tem¬ 
perature  gradients  in  the  presence  of  a  linear  gradient  in  the 
ambient  magnetic  field. 


E.  THESIS  OUTLINE 


3 

Chapter  II  outlines  the  basic  principles  of  the  cryogenic  He 
nuclear  gyro.  The  Bloch  equations,  which  describe  the  dynamics  of  the 
device,  are  presented. 

In  Chapter  III  there  is  given  a  theoretical  analysis  of  cross 
axis  effects.  The  analysis  predicts  the  presence  of  a  false  solution 
when  three  rigidly  mounted  orthogonal  nuclear  gyros  are  used  to  deter¬ 
mine  the  rotational  rate  of  a  body  to  which  the  gyro  triad  is  attached. 
Also  predicted  is  that  certain  oscillatory  rates  can  cause  the  sample 
magnetization  to  align  with  the  applied  magnetic  field.  This  is  impor¬ 
tant  because  this  condition  of  alignment  results  in  no  output  signal. 

An  algorithm  is  developed  to  resolve  the  ambiguity  without  the  need 
for  an  additional  gyro;  and  it  is  shown  that  this  algorithm  may  be  used 
in  conjunction  with  additional  magnetic  field  coils  to  keep  the  sample 
magnetization  orthogonal  to  the  applied  field. 

The  phenomenon  of  kinematic  rectification  is  analyzed  in  Chapter  IV. 
It  is  shown  that,  as  the  coning  frequency  increases,  the  drift 
approaches  a  maximum  value  that  depends  only  on  the  coning  angle.  This 
is  different  from  the  classical  concept  of  kinematic  rectification. 

Chapter  V  describes  an  experiment  to  verify  the  predicted  kinematic 
rectification  effects.  The  salient  features  of  the  experiment  are  that 
it  may  be  conducted  with  a  relatively  simple  apparatus;  and,  by  using 
rotating  magnetic  fields,  no  mechanical  coning  motions  are  required. 

In  Chapter  VI  some  effects  of  magnetic  field  gradients  are  investi¬ 
gated.  Previous  work  [Taber  1978]  had  shown  the  importance  of  magnetic 
field  stability  and  that  relaxation  times  are  reduced  by  the  presence 
of  a  magnetic  field  gradient.  This  analysis  shows  that  a  linear  gradient 
in  the  magnetic  field  also  introduces  a  sensitivity  to  acceleration  and 
temperature  gradients.  These  sensitivities  are  theoretically  predicted 
as  functions  of  sample  cell  size,  temperature,  acceleration,  and  the 
magnitudes  of  the  temperature  and  magnetic  field  gradients. 
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Chapter  II 


CRYOGENIC  He3  NUCLEAR  GYRO  FUNDAMENTALS 

A.  BASIC  PRINCIPLES* 

Consider  a  freely  spinning  body  which  possesses  a  magnetic  dipole 
moment,  p,  Furthermore,  consider  that  p  is  colinear  with  the 
angular  momentum  7  of  the  body  such  that 

d  =  Yj  (2.1) 

where  7  is  a  constant.  Now  place  this  body  in  a  uniform  magnetic 
field  3.  This  will  produce  a  torque  of  p  X  B  on  the  body;  and  the 
classical  equation  of  motion  of  j  is 

I 

j  =  p  X  B  (2.2) 

Or  making  use  of  the  relation  between  3  and  P 

i 

p  =  p  x  -ye  .  (2.3) 

One  might  consider  basing  a  gyro  design  on  such  a  body  in  the  follow¬ 
ing  way.  Place  the  freely  spinning  body  in  a  gyro  housing  and  apply  a 
uniform  magnetic  field  that  is  fixed  in  the  gyro  housing  coordinates. 

Then  if  the  gyro  or  G-coordinate  frame  is  defined  such  that  &  =  B^ib, 
we  can  write  (2.3)  in  the  G-frame  as 


Classical  treatments  of  the  equation  of  motion  of  in  a  magnetic 
field  can  be  found  in  Klttel  [1976]  or  Cohen-T annoud ji  [1977] .  The  same 
results  are  derived  via  a  quantum  mechanical  treatment  [Abragam  1973] . 


G  I 

-  i  -O— I 

n  =  (l  -  u  X  p  .  (2.4) 

If  we  consider  the  case  where  the  gyro  frame  is  nonrotating  with  respect 

•»G— I 

to  an  inertial  frame,  that  is  w  =0,  then  expanding  (2.4)  into 
rectangular  components  in  the  G— frame,  we  get 


d 

dt 


(2.5) 


We  see  that  p  is  uncoupled  from  p  and  p  ,  and  is  constant.  The 
3  A  ^ 

solution  for  p  and  p  is  that  of  a  sinusoidal  oscillator  with  fre- 
A  ^ 

quency  TB3  rad/sec.  A  general  solution  for  p  is  shown  in  Fig.  II-l. 
The  angle,  0,  between  p  and  S  remains  constant  while  p  precesses 
around  B  at  a  radian  frequency  of 


3G 


FIG.  II-l  PRECESSION  OF  A  FREELY  SPINNING  BODY,  POSSESSING  A 
MAGNETIC  DIPOLE  MOMENT  p,  IN  AN  APPLIED  MAGNETIC 
FIELD,  IJ.  The  direction  of  precession  is  that  for 
a  positive  gyromagnetic  ratio  y. 
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Now  for  the  case  where 


-^-1 

CJ 


3G 


we  get,  in  the  G- frame 


’ul" 

0  (>B3+u3)  0  ‘ 

V 

U2 

= 

-(7B3+u)3)  0  0 

U2 

U3 

0  0  0 

U3 

(2.6) 


The  form  of  the  solution  is  the  same  as  (2.5),  but  the  precession  fre¬ 
quency  is  now  (")®3  +  Wg)  rad/sec.  Thus  if  we  know  >8^  and  observe 
either  component  or  p2  with  a  magnetometer,  fixed  to  the  gyro  hous¬ 

ing,  we  can  measure  (>B3  +  ^3)  to  determine  u.'^.  Alternatively,  we 
could  compare  the  phase  of  u  ,  (or  ,  to  that  of  a  stable  oscillator 

with  radian  frequency  7B  and  initially  at  the  same  phase  as  u 

3  1 

(or  Ug).  A  phase  difference  would  then  be  a  measure  of  the  angular 

displacement  of  the  G-frame  with  respect  to  the  I-frame  about  the  3G 

-*G-I 

axis.  The  question  of  how  to  handle  components  of  u  orthogonal 

A 

to  the  3G  axis  is  addressed  in  Ch.  III.  It  is  seen  that  as  long 
as  these  components  are  small  compared  to  7B3,  then  Eq.  (2.6)  is  correct 
to  first  order.  Thus  if  we  can  come  up  with  a  suitable  freely  spining 

body,  a  gyro  housing,  a  magnetometer,  and  a  uniform  magnetic  field,  we 
have  the  potential  for  constructing  a  single  degree-of -freedom  gyro 

that  can  be  used  in  either  a  rate  or  rate-integrating  mode. 

3 

J3.  PROPERTIES  OF  He 

One  such  body  that  posseses  the  requisite  characteristics  is  the 

nucleus  of  the  He  atom.  It  has  intrinsic  spin  a  igular  momentum  of 

4  -1  -1 

l/2  tt,  a  gyromagnetic  ratio  7,  of  -2.038  X  10  rad-eec  -g 

7  1  ^ 23  * 

(-3,244  x  10  Hz-T  ),  and  a  magnetic  dipole  moment  of  1.075  X  10  erg-G 

(1,075  X  10  J-T  ),  Thus  we  cannot  consider  using  only  a  single 

3 

atom  but  rather  a  collection  of  a  large  number  of  He  atoms.  The 
magnetometer  then  senses  the  net  effect  of  the  individual  dipole  moments, 
the  magnetization  of  the  collective  sample,  which  is  defined  as 
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I* 


1 

vol 


(2.7) 


where  N  =  number  of  He  atoms. 

How  then  will  this  collection  of  atoms  behave?  The  individual 
nuclei  cannot  truly  be  considered  freely  rotating  bodies  since  they 
will  influence  each  other  via  dipole  -  dipole  interactions  and  will 
suffer  collisions  with  the  walls  of  the  containment  vessel.  Let  us 
first  look  at  the  equilibrium  condition  of  this  collection  in  a  uniform 
magnetic  field.  Since  the  energy  of  a  dipole  in  an  external  field  is 
-  p  •  if,  then  it  is  energetically  favorable  for  the  dipoles  to  align 
with  the  applied  field.  Thus  we  would  expect  to  find  a  component  of 
IX  along  ^  at  equilibrium.  There  is  no  such  energy  distinction  for 
components  of  orthogonal  to  if  and  thus  at  equilibrium  we  would 

expect  to  have  random  distribution  of  these  components  leading  to  a 
zero  magnetisation  in  that  plane.  Boltzmann  statistics  give  the 
equilibrium  value  of  the  magnetization,  at  temperature  T,  to  be 


%  =  N  |  p|  tanh  .  (2.8) 

Referring  to  (2.4)  and  remembering  that  we  are  interested  in  measuring 
hG-I  a 

the  component  of  w  along  the  3g  axis,  we  see  that  the  equili- 

brium  condition  is  of  no  use  since  p  =  0  here.  So  we  must  consider 

conditions  perturbed  from  equilibrium.  In  particular,  suppose  we 
3 

prepare  the  He  sample  to  have  a  nonzero  magnetization  in  the  1G 
direction.  One  means  of  accomplishing  this  is  via  a  process  called 
optical  pumping  [Colegrove  1963] .  If  we  now  apply  a  field  along  the  & 
axis,  then  those  dipoles  not  aligned  with  the  field  "3  will  process  about 
If  at  the  radian  frequency  yB.  This  is  referred  to  as  Larmor  pre¬ 
cession'  and  yB  is  the  Larmor  frequency.  For  if  uniform,  the  Larmor 
frequency  is  the  same  for  each  of  the  precessing  dipoles;  and  thus, 
if  we  start  out  with  a  nonzero  magnetization,  it  too,  will  precess 
at  the  Larmor  frequency. 
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3 

In  addition  to  the  Larmor  precession,  the  He  sample  will  relax 
back  to  its  equilibrium  state.  If  we  assume  that  the  equilibrium 
state  is  approached  at  a  rate  proportional  to  the  displacement  from 
equilibrium,  then  the  relaxation  of  is  described  by  the  differen¬ 

tial  equation 


d 

dt 


T 

1 


(2.9) 


The  characteristic  time  T  is  referred  to  as  the  longitudinal  relax¬ 
ation  time.  Similarly  for  and  which  have  zero  equilibrium 

values,  we  get 


d 

dt 


(2.10) 


d__ 

dt 


(2.11) 


The  characteristic  time  T2  is  referred  to  as  the  transverse  relaxation 
time,  T  and  T  may  be  nearly  equal;  but,  in  general,  T  is  greater 
than  T  for  reasons  that  will  become  apparent  later. 

Now  combining  both  the  Larmor  precession  and  relaxation  effects,  we 
get  the  Bloch  equations  [Kittel  1976]  in  the  gyro  frame: 


M 

~r  M  =  T0&  X  3),  -  x  Sb,  -  ~  (2.12) 

at  I  l  i  A2 

M 

^M2  =  Y(ti  X  t)2  -  tiP"1  X  M)2  -  ~  (2.13) 

^■M3  =  y(A  X  1)3  -  (S^1  X  M)3  +  .  (2.14) 
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For  a  liquid  sample  of  He  in  He  ,  with  a  solid  hydrogen  wall  coating, 
a  T2  of  140  hours  has  been  obtained.  Thus  for  many  calculations, 
particularly  those  involving  short  term  effects,  the  relaxation  terms 
may  be  neglected.  The  gyro  modeling  performed  in  the  following  chapters 
will  generally  employ  this  assumption. 

£.  MAGNETIC  FIELD  REQUIREMENTS 

3 

One  obvious  difficulty  in  making  a  practical  He  nuclear  gyro  is 

in  obtaining  a  uniform,  static  magnetic  field.  The  magnitude  of  the 

—5 

earth's  magnetic  field  is  on  the  order  of  0.5G  (5  x  10  T)  at  the 

surface  and  varies  in  magnitude  and  direction  from  one  location  to 

3 

another.  So  one  task  is  to  screen  the  He  sample  from  the  ambient 
magnetic  field,  Cabrera  (Cabrera  1975]  has  developed  a  method  of 

_g  —  1 2 

creating  low  field,  10  G  (10  T),  regions  within  a  superconducting 

lead  foil  shield.  This  type  of  shield  has  the  property  of  perfect 

diamagnetism.  That  is,  beyond  a  penetration  depth  which  is  much 

smaller  than  the  foil  thickness,  the  magnetic  flux  density  in  the  foil 

is  zero.  This  is  only  true  provided  the  ambient  field  remains  below 

a  critical  value  and  the  temperature  remains  below  a  critical  value. 

_2 

For  lead,  the  critical  magnetic  field  is  803  G  (8.03  x  10  T)  and  the 
critical  temperature  is  7.2K,  It  is  not  difficult  to  avoid  encounter¬ 
ing  a  critical  field,  but  to  maintain  the  lead  foil  below  7.2K  requires 
that  the  shield  be  kept  in  a  liquid  helium  bath. 

In  addition  to  shielding  from  outside  magnetic  fields,  it  is  also 
required  that  a  high  degree  of  uniformity  of  the  applied  field,  2,  be 
achieved.  If  sixth  order  Helmholtz  coils  are  used  to  generate  3,  then 
gradients  in  will  be  due  mainly  to  the  inhomogeneities  in  the  trapped 

field.  Based  on  current  capabilities  the  major  field  inhomogeneity 

-8  -1 

is  estimated  to  be  a  linear  gradient  of  magnitude  less  than  10  G-cm 
(10-10  T-m_1)  [Taber  1978]. 
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Furthermore  the  field  coils  could  be  made  with  superconducting 
wire  so  that  the  field  will  be  generated  by  the  persisent  supercurrent , 
which  is  inherently  very  stable. 

D.  MAGNETOMETER 
3 

The  remaining  element  of  the  He  nuclear  gyro  is  the  sensor,  the 
magnetometer.  Since  there  is  already  a  requirement  to  maintain  the 
shield  and  the  field  coils  at  cryogenic  temperatures,  a  natural  candi¬ 
date  for  a  magnetometer  is  the  SQUID  (Superconducting  Quantum  Interference 
Device)  magnetometer.  Fortunately,  this  is  the  most  sensitive  magnet¬ 
ometer  currently  available.  Theory  and  use  of  the  SQUID  are  described 

in  such  references  as  Giffard  [1972],  Zimmerman  [1970],  and  Silver  [1967]* 
3 

SQUID  use  for  a  He  nuclear  gyro  is  described  by  Taber  [Taber  1978] . 


* 

See  also  the  S.H.E,  Corp.  Catalog 
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Chapter  III 


CROSS  AXIS  RATE  INPUT  EFFECTS 


A_,  SENSITIVITY  TO  CROSS  AXIS  RATE  INPUT 
3 

Having  a  basis  for  using  a  He  nuclear  gyro  for  measuring  rotations 
about  the  applied  field  axis,  consider  now  the  effects  of  rotations 
about  orthogonal  axes.  Suppose,  as  shown  in  Fig,  III— 1,  that 
B  =  and  w0-*  =  Then  appealing  to  (2.12),  (2.13),  and  (2.14), 

and  neglecting  the  relaxation  terms,  we  get 


it  =  t  +  wI-G  x  S 


(3.1) 


0  7B3  -u2 

VB„  0  0 


(3.2) 


Solving  via  Laplace  transforms,  in  the  gyro  frame 


o  Yb3  -u- 


3(0  «=  Z  (  si  -  -7B  0  0 


“2  0  0 


(3.3) 


(3.4) 


—  sin  ut 
w 


78,  /u2  \2 

- sin  ut  cos  ut  +( — )  (1-cos  ut) 

U  '  U)  ' 


—  sin  ut 
w 


(1-cos  ut) 


(1-cos  ut) 


COS  ut  + 1  11  /  (1-C08  ut) 

'  (ii 


where 


2  ,  „  .2  2 

w  =  (tb3)  +  w2  ■ 


(3.5) 


fttSCEOXiNO  PaCS  bunk-mot  filmed 


FIG.  III-l  VECTOR  DESCRIPTION  FOR  A  GENERAL  CROSS  AXIS  RATE 
INPUT. 


One  effect  of  a  cross  axis  input  is  immediately  apparent.  Even 
with  no  component  of  rotation  about  the  applied  field  (input)  axis, 
the  precession  frequency,  seen  in  the  G-frame,  is  [(yB3)2  +  u>2J  ^  rad/sec. 
If  uncorrected,  the  reading  from  this  gyro  will  be  in  error  by  (w  - 
rad/ sec . 

(w  -  7B3)  =  [(7B3)2  +  w2]  (3.6) 


(3.7) 


Then  expanding  the  term  in  brackets  in  (3.7)  and  keeping  only  the  first 
term  in  (ui^/yfig)  shows  that  the  drift  due  to  a  cross  axis  input 
rate,  is  on  the  order  of  |(u2/t B3)  rad/sec. 
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Figure  III-2  shows  the  effect  of  an  uncorrected  cross  axis  input. 

2 

As  shown  above  the  drift  grows  as  tOg  and  inversely  with  the  Larmor 
frequency,  yBg.  This  suggests  that  cross  coupling  effects  can  be  made 
arbitrarily  small  simply  by  increasing  B  However,  for  a  cross  axis 

7 

input  rate  of  1.0  rad/sec,  the  Larmor  frequency  must  be  10  rad/sec 

—  8 

to  reduce  the  uncorrected  drift  to  5  x  10  rad/sec  (0.01  deg/h). 

3  —2 

For  He  this  requires  that  B  be  about  490  G  (4.9  x  10  T).  The 

O 

gradients  associated  with  such  large  fields  make  them  impractical.  The 
better  approach  is  to  measure  the  cross  axis  inputs  and  then  compensate 
for  their  effects. 

So  let  us  consider  the  general  case  where  there  is  both  an  input 

axis  and  a  cross  axis  component  of  the  input  rate.  Take  the  two  re- 

— >G“  X  a  / 

spective  components  of  w  to  be  w33G  and  m  2G  There  is  no  loss 


There  is  no  loss 


of  generality  here  since  the  direction  of  the  2G  axis,  in  the  plane 
normal  to  the  &  axis,  is  arbitrary.  Proceeding  as  before,  we  get 


°  (^®3+W3)  'W2 


-(  >63+0,3) 


0 


(3.8) 


The  solution  is  thus  the  same  as  (3,4)  with  the  term  783  replaced 

by  (yB  +  u  ).  Thus  the  precession  frequency  seen  in  the  G-frame  is 
J  J  2  2  - 

now  T[(B3  +  w  y/y)  +  (uj^/y)  ]  e’  It  can  also  be  shown  that  the  normal 

to  the  tip  path  plane  (precession  plane)  of  ^  is  parallel  to 

Y(B3  +  u^/y)  &  +  y(0i2/y)  2G.  That  is  the  behavior  which  is  identical 


to  the  case  where  uP-* 


0  and  the  applied  field  is 


B  =  (B  +  -f)3G  +  (-^)  2G.  . 
j  y  y 


This  is  simply  a  statement  of  the  equivalence  between  rotations  and 
magnetic  fields  [Abragam  1973] ,  Thus  the  effect  of  a  cross  axis  rotation, 
<d  ,  is  the  same  as  that  of  apnlying  a  magnetic  field,  oi./y,  in  the  same 


i 


Drift  Rate 
(rad /sec ) 


FIG.  III-2  DRIFT  RATE  VS  LARMOR  FREQUENCY  FOR  A  NUMBER  OF  CROSS 
AXIS  INPUT  RATES,  w  . 

Ct 
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direction  as  the  cross  axis  component  of  <o  .  Thus,  if 

A  A  A 

io^IG  +  +  Wg3G,  the  sensed  precession  frequency,  fi, 

given  by 


-*G-I 

u  = 

will  be 


n 


[(7B. 


2  2  2  1 
+  co„)  +  (0o  +  io.  ] 


(3.9) 


B.  RESOLUTION  WITH  THREE  ORTHOGONAL  GYROS 


As  we  saw  in  general,  the  uncorrected  drift,  [f)-(7B3  +  co3)]  is 
not  tolerable.  So  consider  using  three  identical  gyros,  rigidly  aligned, 
such  that  their  input  axes  are  orthogonal  (Fig.  III-3).  Then  each  gyro 
will  undergo  the  same  rotation  but,  because  the  input  axes  are  orthogonal, 
each  will  be  affected  differently. 

For  the  case  where  to1*-1  =  to, lb  +  io_  £b  +  to  5b,  the  precession 

1  Z  u 

frequencies  sensed  by  each  gyro  are 


It  is  convenient 
2 

(T«3>  ,  giving: 


2 

(7B3 

2  2 

2 

(3.10a) 

“l  " 

rH 

3 

+ 

+  W2 

+  “3 

2 

O 

2 

2 

(3.10b) 

S  = 

+ 

1 

<^3 

+  V 

+  (Og 

2 

2 

2 

<^3  + 

v2 

(3.10c) 

°3  ’ 

wi  + 

W2  + 

(03)  . 

to  normalize  by  dividing  both  sides  of  the  above  by 


jlj  =  (1  +  jo^)2  +  ujg  +  — 3  (3,11a) 

j)2  =  ±>1  +  (l  +  w|)  +  W.3  (3.11b) 

jig  =  jo2  +  +  (1  +  j£g>2  (3.11c) 


where  the  underbars  indicate  normalized  quantities. 


Now  the  question  is,  do  Eqs,  (3.11)  uniquely  determine  w*5-*, 

and,  in  general,  the  answer  is  no.  Knowledge  of  the  time  history  of 
-»b-I 

w  will  not  resolve  the  ambiguity,  as  the  following  example  will 
show.  Suppose  initially  to**-*  =  -l/3  lb  -  l/3  fb  -  l/3  $b.  This  par¬ 
ticular  case  is  resolved  unambiguously  from  (3.11),  Furthermore, 
assume  that  this  condition  persits  for  some  time.  Then  suppose  the 

-»b~  1  AAA 

input  changes  such  that  w  =  -0.3333  lb  -  0.3333  2b  -  0.3333  3b. 


With  only  the  measurements 

are  possible 


solutions  for 


ill,  —2* 


jig  and  knowledge  of  TB^,  two 
One  solution  is  indeed  the  true 


— ^b—  1  A  A 

solution;  the  other,  however,  is  to  =  -0.33337  lb  -  0.33337  2b  - 
0.33337  3b.  The  two  solutions  are  nearly  equal,  though  in  opposite 
directions  from  the  previous  known  solution.  The  smaller  the  step 
away  from  -  l/3  lb  -  l/3  2b  -  l/3  3b  the  closer  are  the  two  solutions. 
Thus  it  is  easy  to  conceive  of  the  possibility  of  locking  on  to  the 
false  solution  from  this  starting  point. 

Let  us  investigate  the  source  of  the  ambiguity  by  considering  the 
space  of  possible  inputs,  w*5”1,  (w  ,  w  ,  w,),  in  which  we  graph  (3-11). 

Zh-f* 

Figure  III-4  shows  the  case  where  to  =0.  Here  the  false  solution  is 


— 1  =  ~2  =  —3  =  1-°* 


We  see  that  the  informa- 


(-2/3,  -2/3,  -2/3),  and 

tion  from  any  one  gyro  defines  a  sphere  on  which  the  possible  values  of 
jo**-1  lie.  The  solution(s)  then  will  be  the  point (s)  common  to  all 
three  spheres.  Figure  III— 4  is  an  isometric  drawing  and  so  the  two 
solutions  plot  as  the  same  point;  nevertheless  one  can  easily  visualize 
the  false  solution  for  this  case.  Note  that  for  three  identical  gyros 
that  the  centers  for  the  three  spheres  are  located  at  (-1,  0,  0), 

(0,  -1,  0),  and  (0,  0,  -1)  and  the  corresponding  radii  are  jl ,  j^, 
and  jig. 

-+b-I 

To  explore  this  situation  further,  consider  the  condition  that  w 
lies  in  the  plane  defined  by  the  centers  of  the  three  spheres.  As  can 
be  seen  from  Fig.  III-5,  there  can  be  no  false  solution.  In  this  plane 
the  solution  is  the  intersection  of  three  circles.  Then,  since  the 
three  centers  are  not  colinear,  there  will  be  at  most  one  point  of 
intersection.  And  since  as  one  moves  out  of  this  plane  in  either 
direction,  the  surfaces  of  the  three  spheres  will  be  moving  away  from 
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one  another,  there  can  be  no  further  points  of  intersection.  So  the 
plane  of  centers  defined  by  (3.12)  is  the  only  region  where 


(3.12) 


the  solution  is  unambiguous. 

What  then  is  the  nature  of  the  ambiguity?  Consider  the  case  where 

ojk-*  lies  outside  the  plane  of  centers.  Choosing  any  two  of  the  spheres 

we  see  that  their  intersection  must  be  a  circle.  Furthermore,  the  plane 

of  centers  cuts  this  circle  in  half;  i.e.,  the  center  of  this  circle 

lies  in  the  plane  of  centers.  Now  since  the  third  sphere  must  also 
-»b~I 

pass  through  u  ,  this  sphere  will  intersect  the  circle  at  two  points. 

One  is  indeed  at  w*5-*  and  the  second  point  is  the  reflection  of  w*5-* 
across  the  plane  of  centers. 

We  can  now  propose  a  simple  criterion  for  resolving  ambiguities. 

X 

Since  one  would  want  to  stay  away  from  the  resonance  condition  (oo  = 
yBg)  in  general,  the  Larmor  frequency  could  be  chosen  such  that 


£  -  T*. 


-*b-I 

for  all  allowable  u> 


Then  the  proper  choice  of  solutions  will  always 


be  on  the  same  side  of  the  plane  of  centers.  Another  possible  criterion, 
though  more  restrictive,  would  be  to  choose 

i*>, i  *  • 


then  the  proper  solution  will  be  the  one  with  the  magnitude  less  than 

i/JS  • 
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£.  SIGNAL  LOSS  DUE  TO  CROSS  AXIS  INPUTS 

The  other  effect  of  a  cross  axis  rotation  is  to  tilt  the  precession 

plane  such  that  it  is  normal  to  the  effective  field,  3  „  „ 

’  ef  f 


ef  f 


“l  ^ 

~  1G  +  —  2G  +  (B3  +  ~)  3G  . 


(3.13) 


The  instantaneous  normal  to  the  precession  plane  will  always  be  in  the 

direction  of  ^  however,  the  manner  in  which  the  normal  tracks 

eft 

depends  on  how  fast  it  changes  direction  and  on  the  initial 

conditions  of  1$  when  the  change  transpires.  If  ^  initially  lies  in 

the  plane  of  precession,  this  condition  can  be  maintained  provided  the 

rate  of  change  of  direction  of  ?  „ .  is  small  compared  to  )B 

eff  '  eff 1 

Violation  of  this  condition  will  result  in  18  tracing  out  a  cone  around 

it  As  an  example,  consider  the  extreme  case  where  initially  It  = 

eff  ’  eff 

=  B^fe;  then  as  M  points  along  the  2G  axis,  there  is  an  instantaneous 
_>G_I  ^ 

change  in  u  to  be  ~u>22G.  This  tilts  the  effective  field  toward  the 

2G  direction.  Now  maintains  a  fixed  angle  with  respect  to  it 

and  thus  moves  up  out  of  the  horizontal  plane  (see  Fig.  III-6).  If  at 

-*  _>G— I 

some  time  later,  when  M  is  at  its  maximum  height,  to  goes  back 

to  zero,  the  precession  plane  will  again  be  a  horizontal  plane;  how¬ 
ever,  no  longer  lies  in  this  plane.  In  fact,  if  the  instantaneous 

-*G_I  -» 

change  in  u  rotated  B  by  the  angle  0,  then  in  the  final 

condition  M  will  make  an  angle  20  with  the  horizontal  plane.  Recall 
that  the  magnetometer  senses  the  horizontal  component  of  M.  Thus  for 
the  example  above  the  signal  has  been  reduced  to  have  a  peak  amplitude 
of  | $ |  cos  20. 

Conceivably  this  pattern  of  events  could  occur  periodically;  and 
would  eventually  become  nearly  coincident  with  the  vertical  axis. 

Such  behavior  was  demonstrated  with  an  analog  simulation  (Fig,  III— 7), 

For  this  simulation  yB  was  chosen  as  1.0  rad/sec.  The  initial 
conditions  were  set  to  have  $(0)  =  |it|  1G  and  u)  =0.  The 
component  of  "S(t)  was  monitored  at  the  analog  console;  then  every 
time  ^  pointed  along  the  £g  axis,  was  made  instantaneously 

Wg^G  (for  this  example  ojg  =  0,05  rad/sec).  Then  when  M  would  point 
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0.0 


FIG.  III-7  ANALOG  SIMULATION  FOR  CROSS  AXIS  INPUTS.  This  is 
a  simulation  of  Eq.  (3.2)  with  yBg  =  1.0  rad/sec 
and  w2  =  0.05  rad/sec  (or  0.0  rad/sec  depending  on 
the  status  of  the  switch).  The  initial  conditions 
are  set  for  $(0)  =  j M {  ffe. 
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~  ->G-I 

along  the  -2G  axis,  w  was  instantaneously  turned  off.  In  this 

manner  $  was  brought  into  alignment  with  the  vertical  axis  in  about 

245  sec  (see  Fig.  III-8).  The  point  here  is  that  oscillatory  behavior 
-»G-I 

of  u)  near  the  resonance  frequency  can  cause  problems.  Remember, 

however,  that  the  second  condition  is  that  the  change  in  direction  of 

^  „  „  be  small  compared  to  ">2  to  avoid  this  problem.  Thus  a  small 

eff  eff  -»G—  I 

amplitude  sinusoidal  oscillation  in  to  at  the  resonance  frequency 

will  not  cause  as  great  a  problem.  Note  also,  though,  that  when  the 

-*G-I 

vertical  component  of  to  is  nearly  equal  and  opposite  to  7B  , 

ytf  ..  can  be  very  small.  Thus  the  resonance  condition  is  not  a 
eff 

constant  but  can  be  reduced  (or  increased)  in  frequency  by  a  Sc  com- 

_>G— I  *  /  1 —  i-»G—  I  | 

ponent  of  u>  .  so  in  addition  to  wanting  yB^  £  l/v3  |max  to 

resolve  ambiguities,  we  have  here  another  reason  to  scale  the  Larmor 

frequency  well  above  the  rotation  rates  we  would  expect  the  device  to 

encounter  in  operation. 


_D ,  USE  OF  ADDITIONAL  MAGNETIC  FIELDS  TO  NULL 
CROSS  AXIS  INPUTS 


It  is  clear  from  Fig.  1 1 1 —  8  that  some  strategy  must  be  employed  to 

keep  the  magnetization  in  the  precession  plane  despite  the  presence  of 

cross  axis  rotations.  We  could  consider  restricting  the  operating  regime 

to  be  always  far  away  from  resonance  conditions;  but  there  may  be  yet 

a  more  satisfactory  approach.  Consider  again  the  Bloch  equations  for 
-*G— I 

an  arbitrary  to 


(7B3+W3)  "W2 


—  M 
dt  2 


-<YB3+u,3) 


(3.14) 


0  M, 


Professor  DeBra  has  suggested  deliberately  introducing  such  a  compo¬ 
nent  by  spinning  the  gyro  about  its  input  axis.  This  increases  the 
effective  Larmor  frequency  without  introducing  the  field  gradient 
problems  that  occur  when  the  applied  field  is  increased. 
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1 1  —  8  MAGNETIZATION  TRAJECTORY  WITH  OSCILLATORY  u>2. 
Note  that  this  simulation  used  a  positive 
7;  the  precession  direction  for  He^  would  be 
opposite  but  the  end  result  the  same.  The 
wiggles  in  the  figure  are  due  to  rough  plotter 
response . 


I f  we  add  coils  for  generating  magnetic  fields  in  the  1G  and  2G  axes 
we  get 


d 

dt 


M 


0 


-(7B3+W3) 


( >B1+w1) 


(TBgtWg) 


M. 


(3.15) 


Thus  if  we  can  provide  two  control  loops  to  keep  both  (7B^  4  )  and 

("yB2  +  u>2>  zero,  we  can  remove  the  effects  of  cross  axis  inputs. 


_E .  THE  THREE  DEGREE-OF-FREEDOM  GYRO 

Since  the  nuclear  gyro  is  sensitive  to  cross  axis  inputs,  it  is 

_*G-I 

theoretically  possible  to  determine  all  three  components  of  w  with 

a  single  device.  To  do  this  three  magnetometers  are  required,  one  to 

-4  -»G—  X 

measure  each  component  of  M.  The  procedure  for  determining  w 

from  measurements  of  the  components  of  ^  is  as  follows: 

1)  Take  three  consecutive  readings  of  iv!  ,  Xf,  ,  and  A?  . 

a  b  c 

2)  Form  the  vectos  aI.  =  and  1ft  =  A?  -]ft 

Aba  B  c  b 

3)  Assuming  uP-1  is  constant  during  the  three  readings,  then 
the  vectors  M  and  M  lie  in  the  precession  plane  and 

>4  — 4  A  B 

Ma  x  Mb/|ma|  lilgl  is  a  unit  vector  in  the  direction 

4  )  1G  -  (")B  4  u^)2G  -  4  Wg)3G  which  is  denoted 

A 

as  n# 


4)  Next  choose  any  of  the  three  measurements,  say  lil  and  deter¬ 


mine  "fl  x 
b 


n  =  \\\ 


sin  t>  or 


sin  4  = 


$  X  n 
b 


PJ 


(3.16) 


5)  The  radius  of  the  precession  plane  is  given  by  [(ft^  |  sin 

\% 


■t  = 


X  n 
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6) 


Now  *|  can  be  determined  (see  Fig.  III-9)  from  the 

radius  of  the  precession  plane,  the  time  between  measurements, 

and  the  measurements  and  si,  . 

a  b 


-1 

a  =  cos 


X  nl2  -  )ttc  - 

2|\  X  n\2 


FIG.  II 1-9  DETERMINATION  OF  w  FROM  THREE  MEASUREMENTS 

OF  lil,  The  view  is  directly  into  the  precession 

Slane,  Note  At  is  the  time  between  measurements 

and  fJ  , 
c  a 
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Assuming  perfect  measurements  and  no  noise  this  scheme  works  very 

well.  Figure  III-10  shows  two  fi?  trajectories  for  a  situation  where 
-*G—  1 

a  step  change  in  w  occurs.  In  this  case  the  size  of  the  step  is 

A 

10%  of  the  Larmor  frequency  and  is  applied  along  the  2G  axis  at  the 
time  when  1  2G.  As  can  be  seen,  the  precession  plane  tips 

down  -  5.7  deg  (0.01  rad)  with  no  control;  but  is  visually  unaffected 
when  a  magnetic  field  is  applied  to  buck  out  the  rotation.  With  control 
applied,  the  vertical  component  of  the  magnetization  reaches  only  a 
value  of  -3.11  x  10  ^  |if(0)  |  ,  and  at  a  sample  rate  of  50  samples  per 
Larmor  period,  the  cross  axis  input  is  nulled  after  only  three  samples. 


fc 


FIG.  III-10  REMOVAL  OF  CROSS  AXIS  EFFECT  WITH  3  DOF  GYRO  SCHEME. 

The  horizontal  precession  plane  is  the  plane  for  zero 
cross  axis  rate,  the  tilted  precession  plane  is  that 
for  u>2  =  0.1  TBgiG.  Use  of  3  DOF  gyro  control  scheme 
removes  effects  leaving  the  initial  horizontal 

precession  plane  virtually  unchanged. 


-32- 


r 


•*w 


There  are  some  problems  with  this  approach.  One  is  that  the 

computation  involved  will  limit  the  sample  rate;  but  probably  more 

crucial  is  the  susceptibility  to  noise.  Figure  III-ll  shows  the  source 

of  the  noise  sensitivity.  Here  the  first  and  third  measurements  are 

-*G-1 

exact  but  the  second  is  noisy.  Since  the  algorithm  to  determine  u 
fits  a  circular  precession  plane  to  these  three  points,  we  see  that  the 
same  amount  of  noise  causes  more  problems  for  fast  sample  rates  than  for 
slow.  So  the  measurement  noise  also  places  restrictions  on  the  sample 
rate. 

One  final  problem  with  this  scheme  should  also  be  mentioned;  and 

that  is  the  effect  of  magnetometer  drift  on  stability  of  the  applied 

field.  Since  the  SQUID  magnetometers  require  a  small,  hut  finite 

current  in  the  pickup  coil,  they  will  also  generate  magnetic  fields.  The 

effects  of  the  pickup  coil  currents,  transverse  to  the  main  field  coil, 

will  be  small;  but  the  pickup  coil  current  in  the  same  axis  will  add 

to  the  applied  field  directly.  If  this  current  were  constant,  its  effect 

-9 

could  be  modeled  out;  however,  a  current  change  of  lxio  amps  gives  a 
drift  of  ~5xio-5  rad/sec  for  the  device  in  Ch.  V. 

F.  THREE  ORTHOGONAL  SINGLE  DEGREE-OF-FREEDOM  (SDOF)  GYROS 

We  now  return  to  see  what  we  can  do  with  three  orthog  1  SDOF  gyros 
each  having  one  or  two  magnetometers  orthogonal  to  each  other  (if  two) 
and  to  the  applied  field  axis.  Since  we  now  have  added  field  coils  in 
the  cross  axes  to  null  cross  axis  rates,  this  is  not  the  same  problem 
as  was  investigated  in  Sect.  Ill— B,  These  added  field  coils  take  out 
different  components  of  the  body  rate  for  each  of  the  three  gyros.  So 
while  the  solution  for  w*5-1  is  still  the  Intersection  of  three  spheres 
the  centers  of  the  three  spheres  are  no  longer  fixed.  Just  as  the 
applied  field  moves  the  center  of  a  given  sphere  along  the  negative 
input  axis  for  that  particular  gyro,  so  too  do  the  feedback  fields  move 
the  center  along  the  cross  axes. 

The  question  of  ambiguities  is  not  as  clear  here  as  it  was  for  the 
case  with  no  cross— axis  field  coils.  However,  one  simple  algorithm 
appears  to  work  very  well  for  determining  w  from  the  three  gyros. 


-33- 


NOISY  MEASUREMENT 


ACTUAL  PRECESSION 
PLANE 


FIG.  III-ll  EFFECT  OF  MEASUREMENT  NOISE  ON  3  DOF  GYRO  SCHEME. 

The  upper  sketch  shows  the  effect  of  a  noisy  meas¬ 
urement  for  a  sample  rate  with  respect  to  precession 
rate.  The  computed  precession  plane  is  orthogonal 
to  the  actual  precession  plane.  The  lower  sketch 
shows  the  effect  of  the  same  amount  of  noise  when  the 
sample  rate  is  slow.  Here  the  actual  and  computed 
precession  planes  are  nearly  the  same. 


The  procedure  is  as  follows. 

1)  Measure  the  precession  frequency  for  each  gyro. 

2)  Compare  this  measured  frequency  with  the  zero  input  Larmor 
frequency  and  take  the  difference  to  be  the  rate  about  the 
input  axis. 

3)  Use  these  estimated  rates  to  determine  the  required  cross 
axis  fields. 

4)  Continue  this  procedure. 

Figure  III-12  shows  the  response  to  a  step  in  To'3-*,  The  step 
size  is  25%  of  the  Larmor  frequency  and  is  applied  along  the  3b  axis. 
The  significance  to  the  time  scale  on  Fig.  III-12  is  that  the  preces¬ 
sion  frequency  of  each  gyro  is  assumed  to  be  determined  in  one  time 
unit.  So  let  us  look  at  what  is  involved  in  determining  the  precession 
frequency  from  the  SQUID  magnetometer (s)  output.  If  one  assumes  that, 
during  the  measurement  period,  the  precession  frequency  is  constant  and 
the  gyro  housing  rotation  rates  are  small  compared  to  the  Larmor  fre¬ 
quency,  then  the  SQUID  output  will  be  a  sinusoidal  signal,  within  a 
certain  band,  centered  at  the  Larmor  frequency.  If  the  unwanted  low- 
and  high-frequency  bands  are  filtered  out  then  the  SQUID  output  Mg(t) 
will  have  the  form 

M(t)  =  ^  sin  wt  +  B  cos  wt  .  (3,17) 

s 

Consider  sampling  this  signal  at  a  fixed  sample  period,  T  sec.  Then 
it  is  shown  that  the  nth  sample  is  given  by 

M  (n)  =  B6(T)  +  (A  sin  uT  -  B  cos  u>T)6(T-l) 

s 

+  2  cos  u ft  Mg(n-1)  -  Mg(n-2)  .  (3.18) 

We  can  write  this  out  in  matrix  form  as 
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f— 

1 

M  (0) 

i 

0 

I 

0 

0 

s 

1 

Ms(1) 

0 

1 

1 

M  (0) 

0 

— 

— 

1 

— 

—  _  - 

B 

M  (2) 

0 

0 

1 

M  (1) 

M  (0) 

s 

= 

1 

s 

S 

A  sin  (oT  - 

B  cos  wT 

M  (3) 

0 

0 

1 

M  (2) 

M  (1) 

S 

1 

1 

s 

s 

2  cos 

toT 

M  (4) 

0 

0 

1 

M  (3) 

M  (2) 

s 

1 

s 

s 

-1 

• 

• 

• 

1 

• 

L  _ 

1 

The  upper  left  partition  is  due  to  initial  conditions  so  if  we  wait 
two  sample  periods  we  have 


M  (2) 
s 

0— 

Ms(l) 

M  (0) 
s 

Ms(3) 

II 

Ms(2) 

Ms(l) 

Ms(4) 

Ms(3) 

Ms(2) 

w  ■  - 

V _ _ 

J 

y  a 


2  cos  wT 

-1 


(3.20) 


Now  define  2  cos  uT  =  a,  and  -1  =  a„  and  consider  the  case  where  we 

1  2 

have  noise  v  such  that 


y  =  Aa  +  v  .  (3.21) 

We  know,  via  measurements,  both  y  and  A;  and  we  can  get  a  best  least 
squares  estimate  of  a  as 

a  =  (ATA)_1ATy  .  (3.22) 

Since  we  know  a_  =  -1  we  could  consider  the  convergence  of  a  to  -1 

a  4 
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as  a  simple  test  for  a  good  estimate  of  a^.  Finally,  we  estimate  the 
precession  frequency  from  as 


A 

D 


Then  the  estimates  of  gj,  one  for  each  gyro,  can  be  used  to  determine 
the  fields  necessary  to  cancel  the  cross  axis  rates.  Note  that  if 
fields  are  generated,  there  will  be  changes  in  the  precession  frequencies 
of  the  three  gyros.  Thus  the  next  two  samples  will  again  exhibit 
the  transient  terms  in  (3,19). 
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Chapter  IV 


KINEMATIC  RECTIFICATION 


A.  STATEMENT  OF  THE  PROBLEM 

One  rectification  mechanism  has  been  discussed  in  Ch.  Ill,  i.e., 
the  rectification  of  cross  axis  rates.  It  was  seen  that  a  cross  axis 
rate  increased  the  precession  frequency  independent  of  the  sign  of  the 
cross  axis  rate.  Thus  oscillatory  cross  axis  rate  inputs  will  produce 
an  average  drift  even  if  the  average  cross  axis  rate  is  zero.  Another 
well  known  rectification  drift  for  SDOF  gyros  is  kinematic  recti¬ 
fication  [Cannon  1957"],  The  mechanism  for  kinematic  induced  drift, 

3 

as  it  applies  to  the  He  nuclear  gyro  is  investigated  here. 

Consider  the  situation  shown  in  Fig.  IV-1.  The  I-frame  is  an 

3 

inertial  frame  and  the  G-frame  is  the  He"  gyro  frame.  Assume  that  the 
applied  field  is  along  the  3G  axis  and  that  0^  and  0^  are  driven 
such  as  to  maintain  a  constant  angle  [?  between  the  3G  and  31  axes 
as  the  3g  axis  cones  around  the  31  axis  at  a  frequency  rad/sec. 

Then 

_*G- 1  *  a  •  r 

^  =  ?1  1G  +  02  21  (4.1) 

wG_I  =  ^10  +  f?2  cos  0  £g  -  02  sin  ^  3G  .  (4.2) 

Thus  if  the  gimbal  angles  are  given  by 

0_(t)  =  -P  cos  w  t  (4.3a) 

1  z 

0_(t)  =  P  sin  w  t  (4.3b) 

Ci  Z 

then  for  small  angles  sin  0^  •»=  0^  and  we  see  the  component  about  the 
3g  axis  is  approximately  given  by 

•  *22 

-0  sin  0  ea  -0  0  =  -p  10  cos  w  t  .  (4.4) 

Z  1  ci  x.  Z  Z 
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INPUT  AXIS  CONING  OF  A  SINGLE  DEGREE  OF  FREEDOM  (SDOF) 
GYRO.  Angles  0  and  9  are  driven  to  cause  the 

*■  ^  A 

axis  to  cone  about  the  31  axis,  at  a  frequency  ioz , 
while  maintaining  a  constant  angle  g  between  the  two 
axes . 


This  component  has  an  average  value  despite  the  fact  that  the  gyro 
has  no  average  displacement  about  the  3G  axis  with  respect  to  the  iner¬ 
tial  frame.  Since  the  gyro  is  intended  to  measure  the  ic  component  of 
-»G— I 

co  though  we  expect  to  see  (-0-0  )  as  the  gyro  output.  This  is 

^  L  3  Vg 

the  drift  due  to  kinematic  rectification  and  is  given  by 


(-eo0. ) 

2  1  avg 


-  \  d2 

peak  peak 


1  2 

=  '2^2 


(4.5) 

(4.6) 


And  so  this  term  is  seen  to  depend  on  the  cone  half  angle  p  and  the 

frequency  of  the  coning,  u;  .  Note  that  nothing  has  been  assumed  about 

z 

the  nature  of  the  gyro  itself  other  than  it  is  a  single  degree-of-f reedom 
gyro.  Thus  one  would  expect  to  see  this  behavior  in  all  single  DOF  gyros. 
Indeed,  it  has  been  observed  in  conventional  mechanical  gyros  and  also 
laser  gyros.  To  date  no  attempt  has  been  made  to  observe  kinematic 
rectification  in  nuclear  gyros.  Chapter  V  will  describe  an  experiment 
design  to  make  this  observation. 

_B.  ANALYSIS  OF  KINEMATIC  RECTIFICATION  EFFECTS  FOR  THE 

He3  GYRO 

We  can  solve  for  the  effects  of  coning  motion  for  the  nuclear  gyro 
using  the  Bloch  equations  and  assuming  relaxation  effects  can  be  neg¬ 
lected.  This  problem  is  most  easily  solved  in  a  rotating  frame  where  the 
applied  field,  "8  =  B  5g,  is  constant  (see  Fig.  IV-2).  This  rotating 

«J 

/■' 

or  R-frame  rotates  at  the  coning  frequency,  u>z,  about  the  3R  axis. 

So  in  this  frame  the  applied  field  is  constant  and  is  given  by 

iL  =  B  sin  ft  1R  +  B  cos  p,  $R  (4.7) 

K  j  o 

=  B  fR  +  B  §R  .  (4.8) 

1  z 


-41- 


Then  the  equation  of  motion  for  M  in  the  R-frame  is 


R  I 

«l  =  it  +  UI-R  x  ft  . 

Writing  out  the  components  in  the  R-frame  we  have 


’Ml‘ 

0 

(7B  -u>  ) 
z  z 

0 

M 

1 

d 

dt 

M2 

= 

1 

5 

n, 

nE 

0 

*1 

U2 

M3 

L  0 

-*1 

0 

-M3. 

(4.9) 


(4.10) 


Finally  we  solve  for 


M  in  the  R-frame  to  get 


M(t) 


(4.11) 


3(t)  = 


/>*  i\2 

cos  wt  +  (1-cos  ut)  - - —  sin  wt 


u 


u 

z  z 


7B.  (7B  -u  ) 
1  z  z 


sin  ut 


(1-cos  ut) 


COS  ut 


78, 


YB,(7B  -u  ) 

1  z  z 


(1-cos  ut) 


7B, 


sin  ut 


- -  sin  ut  cos  ut  +( - — }  (1-cos  ut) 

u  \  u  / 

(4.12) 


>1(0) 


where 


w  = 


(7fix) 


<*»«  - 


"*) 


(4.13) 


Now  if  Fig.  IV-2  depicts  the  initial  conditions,  then  $(0)  =  -  |$|  2h  and 
the  solution  is 
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i 


OB  -  w  ) 

M^(t)  =  -  |^|  - — - —  sin  cot 

co 


M  (t)  =  -  j  M |  cos  cot 

M  (t)  =  111  I  — —  sin  cot  . 

3  1  1  to 


We  now  make  use  of  the  transformation 


/V 

“ 

- 

11 

sin  co  t  -cos  to  t  0 

IE 

z  z 

A 

21 

_ 

coscot  sin  to  t  0 

£r 

z  z 

/N 

/N 

31 

0  0  1 

3R 

-  . 

- 

- 

(4.14a) 


(4.14b) 

(4 . 14c  ) 


(4.15) 


the  solution  in  I-frame  coordinates  becomes 


Mx(t)  =  -  |M 


M2(t) 


L  w 
r  ob -w ) 

-  I*l[ - „ -  sin 


sin  cot  sin  co  t  -  cos  cot  cos  co  t 
z  z 


cot  cos  co  t  +  cos  cot  sin  co  t 
z  z 


_  TB1 

M^(t  )  =  |  M |  — —  sin  cot 


(4.16a) 


(4.16b) 


(4.16c) 


Now  if  we  look  at  the  M^(t)  component  in  the  I-frame  we  can  simplify 
(4.16a)  to  get 


v0 = - 


)B  -co  -co  ya  -co  +co 

2  2  Z  Z 

-  COS(U)-OJ  )t  -  —  1  cos(u)  + 

(0  Z  CO 


U)z)t}  • 


(4.17) 


In  as  much  as  the  1G  axis  never  gets  far  from  the  il  axis,  then 
if  the  coning  motion  produced  no  effect,  a  magnetometer  in  the  fG  axis 
would  see  a  signal  that,  on  the  average,  looked  like  jl8 1  cos  yB^t. 
Equation  (4.17),  however,  shows  that  this  is  not  the  case.  The  question 
now  is  does  the  deviation  from  | ft! j  cos  yB^t  agree  with  drift  term  we 
would  expect  from  kinematic  rectification.  Table  IV-1  shows  a  variety 
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of  conditions  comparing  the  drift  rate  calculated  from  (4.17)  with  that 
predicted  from  (4,6). 


Table  IV-1 

COMPARISON  OF  DRIFT  PREDICTED  BY  ANALYSIS  WITH 
CLASSICAL  KINEMATIC  RECTIFICATION  DRIFT 


TBg  =  1.0  rad/sec; 

p  =  0.5  deg 

U) 

z 

Drift  Predicted  by 

Eq.  (4.17) 

Drift  Predicted  by  Classical 

Kinematic  Rectification  (Eq.  4.6) 

(rad/sec) 

(rad/sec) 

(rad/sec ) 

0 

0 

0 

-4 

-9 

-9 

1  a  10 

4.0  x  10 

4.0  x  10 

,  -3 

-8 

-8 

1  X  10 

3.8  X  10 

3.8  X  10 

,  -2 

-7 

-7 

1  A  10 

3.85  x  10 

3. 81  X  10 

-2 

-6 

-6 

5  x  10 

2.00  X  10 

1.90  X  10 

-1 

-6 

-6 

1  X  10 

4.23  X  10 

3.81  X  10 

-1 

-5 

-5 

5  a  10 

3.81  X  10 

1.90  X  10 

-1 

-5 

-5 

o 

< 

8.88  X  10 

2.67  X  10 

-i 

-4 

-5 

9  a  10 

3.42  X  10 

3.43  X  10 

-1 

-3 

-5 

9.9  a  10 

3.24  x  10 

3.77  x  10 

TB3  cos  P 

-3 

8.69  x  10 

-5 

3.81  X  10 

l 

_ _ _ _ 1 

Table  IV-1  shows  good  agreement  between  the  analysis  and  the  classical 

kinematic  rectification  drift  up  to  an  w  which  is  about  10%  of  the 

z 

Larmor  frequency.  For  wz  >  O.IYB^,  the  drift  predicted  from  (4.17) 
is  greater  than  would  be  predicted  from  the  kinematic  rectification 
effect  alone.  Part  of  this  increase  is  due  to  cross  axis  rate  rectifi¬ 
cation.  Since  the  precession  rate  is  the  square  root  of  the  sum  of  the 
squares  of  each  component,  we  must  consider  the  average  squared  values 
of  the  cross  axis  rates.  So  combining  these  two  rectification  mechanisms, 
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we  could  predict  the  magnitude  of  the  precession  rate,  PR,  to  be 


PR  =  y  (>B3  +  if2./  r  |(f./  +  |(fto/  .  (4.18) 

By  including  the  cross  axis  rate  rectification  effect,  we  do  get  better 
agreement  with  the  analysis;  but  even  in  this  case  agreement  is  not 
good  as  we  approach  the  resonance  condition  =  "yB^  cos  p)  .  Table 

IV-2  summarizes  these  results. 


Table  IV-2 

COMPARISON  OF  DRIFT  PREDICTED  BY 
EQS .  (4.17)  and  (4.18) 


■yBg  =  1.0  rad/sec; 

f  =  0.5  deg 

Wz 

Drift  from  (4.17) 

Drift  from  (4.18) 

(rad/sec ) 

(rad/sec ) 

(rad/sec) 

-4 

-9 

-9 

1  X  10 

4.0  X  10 

4.0  X  10 

-3 

-8 

1  X  10 

3.8  X  10 

3.8  X  10 

-2 

-7 

-7 

i  x  io 

3.85  X  10 

3.85  X  10 

-l 

-6 

-6 

1  X  10 

4.23  X  10 

4.19  X  10 

-1 

-5 

-5 

5  x  10 

3,81  X  10 

2.86  X  10 

To  understand  why  the  drift  rate  given  by  (4.17)  is  still  greater 
than  that  predicted  by  (4.18)  as  resonance  is  approached,  consider  what 
happens  at  resonance.  At  the  resonance  condition  the  rotating  component 
of  the  applied  field  rotates  about  the  Si  axis  at  the  same  rate  that 
the  magnetization  vector  precesses  about  it.  For  the  initial  conditions 
shown  in  Fig.  IV-2,  13  will  precess  around  the  stationary  component  of 

the  field  at  a  rate  TB^  cos  @  rad/sec.  Since  also  remains  perpen¬ 
dicular  to  as  13  precesses  about  the  static  field,  will  also 
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precess  about  B^  at  a  rate  YBg  sin  3  rad/sec.  These  two  precessions 
occur  simultaneously  so  both  sum  and  difference  frequencies  are  observed. 
Coming  into  resonance  for  the  particular  initial  conditions  chosen 
here  we  see  the  sum  of  the  two  frequencies,  namely  cos  3  +  yB^  sin  3) 

rad/sec,  which,  for  the  conditions  used  to  generate  Tables  IV-1  and  IV-2, 
is  1.00869  rad/sec.  Thus  we  get  the  drift  term  8.69  x  10  rad/sec. 

Had  we  chosen  for  our  initial  conditions  a  situation  which  had  and  13 

parallel,  we  would  have  a  precession  frequency  of  yB^  cos  3  rad/sec  at 
resonance,  because  the  rotating  B^  would  have  no  influence  on 

£.  FREQUENCY  DEPENDENCE 

At  this  point  we  have  built  a  good  case  for  a  classical  kinematic 

3 

rectification  drift  effect  for  the  He  nuclear  gyro.  We  have,  however, 

only  examined  cases  for  which  was  in  the  same  direction  as,  and  no 

greater  than,  the  Larmor  precession.  As  we  pass  beyond  the  resonance 

condition  another  mechanism  takes  over.  Referring  to  (4.17),  it  may  be 

shown  that  from  u  =  0  to  just  below  resonance,  the  coefficient  of 
z 

the  cos(io  +  w  )t  term  is  approximately  2  while  that  of  the  cos(w-w  )t 
z  z 

term  is  nearly  zero.  The  opposite  is  true  beyond  resonance  and  the 
former  is  nearly  zero  while  the  latter  is  about  -2  (at  resonance  the 
first  is  -1.0;  the  second  is  1.0).  Table  IV-3  shows  the  drift  for 
coning  frequencies  above  the  resonance  condition.  Note  that  coming  out 
of  resonance  (u^  =  1.0  rad/sec)  the  precession  frequency  is  nearly 
(yB3  cos  3  -  yB3  sin  3)  which  would  give  a  drift  rate  of  -8.76  x  10  3 
rad/sec. 

The  region  above  resonance  is  characterized  by  a  drift  rate  opposite 

in  direction  to  the  drift  rate  below  resonance  and  asymptotically  approach- 

—5 

ing  -3.81  x  10  rad/sec.  It  may  seem  strange  that  the  drift  rate  should 
change  signs  when  both  the  kinematic  rectification  and  cross  axis  rate 
rectification  effects  predict  a  steadily  increasing  positive  drift  rate 
with  increasing  However,  what  has  happened  is  that  these  effects 

no  longer  apply.  Recall  from  the  discussion  in  the  previous  chapter  that 
when  the  direction  of  the  effective  field  changes  slowly,  the  precession 
plane  is  able  to  follow,  and  remains  in  the  precession  plane. 
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Table  IV-3 


f 


DRIFT  RATES  FOR  to  ABOVE  RESONANCE 
z 


YB3  =  1.0; 

3  = 

0.5  deg 

D 

Z 

Drift 

Rate 

From  (4.17) 

(rad/sec ) 

(rad/sec ) 

1.0 

-8.73 

X 

IQ'3 

5.0 

-4.76 

X 

10-5 

10.0 

-4.23 

X 

10~5 

50.0 

-3.  89 

X 

10-5 

100.0 

-3.85 

X 

10-5 

500.0 

-3.82 

X 

10-5 

3 

When  this  is  the  case  the  He  nuclei  behave  as  if  they  were  gimballed, 

and  their  spin  axes  follow  the  applied  field  as  it  cones  about.  In 

this  regime  (4,18)  applies.  As  the  coning  frequency  increases,  however, 

the  magnetization  vector  cannot  stay  in  the  precession  plane.  Then  at 

3 

sufficiently  high  frequencies  the  He  nuclei  do  not  respond  to  the  rotat¬ 
ing  component  of  the  applied  field  at  all,  and  they  become,  in  effect, 
uncoupled  from  the  coning  motion  of  the  gyro  case.  Why  then  the  drift? 

The  answer  is  that  since  the  applied  field  is  made  up  of  both  a  rotating 

3 

and  a  static  component;  and  for  uz  »  7B3  cos  p  the  He  responds  only 
to  the  static  component,  the  precession  frequency  will  approach  cos  p. 

Thus  it  is  slower  than  the  Larmor  frequency,  and  in  the  limit  approaches 

a  drift  rate  of  (7B„  cos  P  -  yB  ) ,  which  for  the  conditions  in  the  tables 

-5  J ,  3 

is  -3.81  x  10  rad/sec. 

Figure  IV-3  shows  the  magnitude  of  the  drift  vs  the  coning  frequency 
for  =  1,0  rad/sec  and  P  =  0.5  deg.  The  solid  curve  is  for 

opposite  in  direction  to  the  Larmor  precession,  and  the  dashed  curve  repre¬ 
sents  the  case  where  the  directions  are  the  same.  Note  that  the  resonance 
condition  applies  only  in  the  latter  case, 
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DRIFT  VS  CONING  FREQUENCY 


We  can  see  at  the  low  coning  frequencies,  s  0.5  |tB3|,  that 

the  dominant  effects  are  the  two  rectification  mechanisms;  Eq.  (4.18) 
gives  good  results  in  this  range.  Then  as  increases  in  magnitude 

either  resonance  or  uncoupling  from  the  motion  sets  in  and  (4.17)  must 

3 

be  used.  Finally,  in  the  limit  of  large  the  He  nuclei  do  not 

respond  to  the  B  component  of  the  field  and  the  drift  rate  approaches 
—5  ^ 

-3.81  x  10  rad/sec.  Figure  IV-4  summarizes  the  regions  of  influence 
for  the  case  where  (3  =  0.5  deg. 
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Chapter  V 


AN  EXPERIMENT  TO  MEASURE  KINEMATIC 
RECTIFICATION  EFFECTS 


A-  objective 

Because  of  the  nonconventional  behavior  with  regard  to  kinematic 
rectification,  it  was  desired  that  some  experimental  verification  be 
performed.  In  particular,  the  experiment  should  be  designed  to 
accomplish  the  following, 

a)  Verify  the  frequency  dependence  of  both  kinematic  and  cross 
axis  rate  rectification  mechanisms  by  varying  the  coning  fre- 
quenc  y. 

b)  Separate  kinematic  from  cross  axis  rate  rectification  effects 
by  changing  the  sign  of  u  . 

c)  Verify  the  uncoupling  from  the  coning  motion  at  high  w^. 

_B.  LABORATORY  FACILITIES 

As  discussed  previously,  Taber  [Taber  1978]  had  performed  extensive 

3 

experiments  measuring  the  longitudinal  relaxation  times  of  polarized  He  . 

3 

Thus  there  already  existed  a  cryogenic  facility  for  performing  He 
nuclear  gyro  experiments  (see  Fig,  V-l).  Furthermore,  the  cryostat 
probe  used  for  the  relaxation  time  experiments  (Fig.  V-2)  was  still 
available , 

Longitudinal  relaxation  times,  T  of  140  hours  had  been  obtained 

1  3 

with  this  device  [Taber  1978]  using  a  liquid  sample  of  0.07%  He  in 
4 

He  solution.  Transverse  relaxation  times,  T  ,  have  not  been  measured; 
but  since  we  will  now  have  a  lower  ambient  magnetic  field  (~  2  x  10  G 
or  2  x  10  ■12t)  and  will  remove  the  magnetic  contaminant  that  was  present 
for  the  earlier  T  measurements,  we  can  expect  T  to  be  of  the  same 
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MOUNT  FOR  OPTICAL 
PUMPING  &  GAS 
HANDLING  APPARATUS 


FIG.  V-l  GENERAL  EXPERIMENTAL  CONFIGURATION.  Not  shown  are  the 
field  coils  for  the  optical  pumping,  the  optical  pump¬ 
ing  apparatus  and  the  gas-handling  apparatus.  [From 
Taber  1978.] 
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5cm 


FIG.  V-2  SCHEMATIC  DEPICTION  OF  THE  BOTTOM  PORTION  OF  THE 
CRYOSTAT  PROBE  [from  Taber  1978]. 


-55- 


order  of  magnitude  as  the  T  measurement.  Further,  by  carefully 

1  3  4 

metering  the  flow  of  the  polarized  He  -  He  mixture  into  the  sample 
cell,  a  uniform  cell  wall  coating  of  frozen  hydrogen  can  be  maintained. 
This  will  reduce  the  relaxation  due  to  wall  collisions  [Taber  1978] 
thus  increasing  both  T  and  So  it  looks  very  promising  that  T ^ 

should  be  sufficiently  long  to  perform  the  desired  experiments. 

There  were,  however,  some  problems  with  the  existing  experimental 
configuration.  The  first  was  that  the  probe  had  only  two  orthogonal 
field  coils.  The  second  was  that  the  material  used  to  construct 
the  coil  forms,  delrin,  was  not  sufficiently  free  of  ferromagnetic 
contaminants.  Both  of  these  problems  were  solved  by  redesigning  that 
portion  of  the  device.  Corning  Macor,  a  machinable  glass  ceramic,  was 
chosen  as  the  material  for  the  coil  forms.  Samples  of  this  material 
were  found  to  be  magnetically  cleaner  than  similar  samples  of  delrin 
[Cabrera  1975], 

One  final  problem  was  that  for  the  given  experimental  configuration 
mechanical  inputs  are  difficult.  The  gyro  is  located  in  the  cryostat 
probe  which  is  lowered  into  the  dewar;  and  only  limited  rotation 
around  the  vertical  axis  is  possible. 


£.  BASIC  CONCEPT 

Since  a  controlled  mechanical  coning  motion  is  not  feasible,  the 
experiment  can  be  performed  utilizing  the  equivalence  between  mechanical 
rotations  and  applied  magnetic  fields.  Thus  instead  of  physically 
coning  the  gyro  input  axis,  the  rotating  component  of  the  magnetic 
field  is  generated  with  the  coils  in  the  fG  and  2G  axes.  Then  since 
the  gyro  will  not  be  subject  to  any  rotations  (except  earth  rotation) 
we  can  consider  it  to  represent  the  I-frame;  and  thus  a  magnetometer 
located  in  the  fb  axis  will  sense  the  M^(t)  defined  in  (4.16a). 

Now  if  we  express  I?  in  the  I-frame  we  get 

B*  =  | B  |  sin  p  sin  w  t  fl  +  ] B  ]  sin  p  cos  w  t  21 

Z  J  z 

(5.1) 


+  | B  |  cos  p  §1  . 


1 
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Thus  for  a  given  cone  half  angle  (3,  a  coning  frequency  to  and  a  field 
magnitude  B^,  we  can  generate  ^  with  three  orthogonal  coils. 

For  the  purpose  of  this  experiment  the  So-axis  or  static  field  coil 

—  6  —  1 0 

was  chosen  to  be  50  X  10  G  (50  X  10  T),  This  will  give  a  Larmor 

precession  of  about  1  rad/sec.  Two  factors  which  influenced  this 
choice  were  first  the  field  gradients  go  up  with  the  field  magnitude 
and  second  such  a  slow  precession  frequency  allows  for  switching  of 
various  fields  manually. 

Next,  in  order  to  make  the  drifts  due  to  the  coning  motion  dominant 
(and  large  enough  to  see  on  a  strip  chart  recording),  p  was  chosen 
large  enough  to  give  an  effect  about  an  order  of  magnitude  greater  than 
earth  rate.  So  choosing  g  =  20.35  deg  and  =  1.0  rad/sec  gives 

=  0.371  rad/sec.  Now  Table  V-l  shows  a  comparison  of  the  calculated 
drift  (Eq.  4.17)  rate  and  that  predicted  by  the  classical  kinematic 
rectification  effect.  All  of  the  drift  rates  shown  in  Table  V-l  should 
be  large  compared  to  any  other  uncalibrated  drift  rates.  They  are  also 
large  enough  to  observe  in  a  time  that  will  be  small  compared  to  the 
expected  relaxation  time. 


D.  FIELD  COIL  DESIGN 


Since  magnetic  field  uniformity  is  critical,  all  four  coils  (includes 
the  magnetometer  pickup  coil)  are  designed  in  a  helmholtz  configuration 
(see  Fig.  V-3).  For  this  configuration,  where  each  coil  has  N  turns 
and  the  same  current,  I,  the  total  magnetic  field  at  any  point  on  the 
z— axis  is  given  by 


B  (z)  = 

z 


^Nlr 

1 

2 

71  2.3/2 

_(r  +z  )  ' 

(5.2) 


where  y,0  is  the  permeability  of  free  space. 
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Table  V-l 


COMPARISON  OF  CALCULATE!)  DRIFT  RATES  WITH  THOSE 
PREDICTED  BY  CLASSICAL  KINEMATIC  RECTIFICATION  EFFECTS  ALONE 


*3 

=  1.067  rad/sec  p  = 

20.35  deg 

Wz 

Calculated  Drift 

Kinematic  Rectifica— 

Eq.  (4.17) 

tion  Drift  Eq .  (4.6) 

(rad/sec) 

(rad/sec) 

(rad/sec ) 

0.  01 

6.299  X  10-4 

-4 

6.037  x  10 

0.05 

3.269  X  10-3 

3. 154  x  10“ 3 

0. 10 

6.863  X  10-3 

6.307  x  10“ 3 

0.50 

5.596  X  10-2 

3.154  x  10" 2 

0.  70 

1.104  X  10-1 

4.415  x  10“2 

0.90 

2.174  x  10_1 

5.677  x  10“2 

-  0.01 

-6.185  x  10~4 

-6.307  x  10~4 

-  0.05 

-2.985  X  10-3 

-3.154  x  10“3 

-  0. 10 

-5.721  x  10-3 

-6.307  x  10-3 

-  0.50 

-2.140  x  10-2 

-2  i 

-3.154  x  10 

-  1.00 

-2 

-3.247  x  10 

-6.307  x  10“ 2 

-  5.00 

-5.513  x  10-2 

-3.154  x  10_1 

-10.0 

-6.034  x  10-2 

-6.307  x  10-1 

-50.0 

-6.525  x  10-2 

-3. 154 

—4 

Note:  The  smallest  expected  drift  is  -6.185  x  10  rad/sec  = 

-2.126  deg/minute. 


It  is  easy  to  show  from  (5.2)  that  both  the  first  and  second  deriv- 
itives  of  the  field,  with  respect  to  z  ,  vanish  at  the  center  (z  =  r/2). 

The  values  for  "YB^  and  were  given  in  the  previous  section  as 

1.0  rad/sec  and  0.371  rad/sec  respectively.  Thus  the  3G  axis  coil 

“*  6  “10 

must  produce  a  steady  50  x  10  G  (50  X  10  T)  field,  while  the 
1G  and  2G  axis  coils  must  generate  peak  fields  of  18.199  x  10  6G 
(18.199  x  10  10T).  To  accomplish  this,  these  coils  were  designed  as 
specified  in  Table  V-2.  Figure  V-4  shows  the  calculated  field  profile 
across  the  sample  cell  for  each  coil  pair.  As  mentioned  earlier  the  coil 
forms  are  machined  out  of  Macor,  After  machining  the  Macor,  forms  are 
degreased  and  processed  in  an  acid  bath  to  remove  any  possible  magnetic 
contaminants  picked  up  in  the  machining  process.  Assembly  is  performed 
on  a  clean  bench.  Figure  V-5  shows  the  coil  form  assembly. 

Table  V-2 

FIELD  COIL  SPECIFICATIONS 


Gyro  Axis 

Diameter 
(meters ) 

Number  of  Turns 

Peak  Current 
(amps ) 

1G 

4.0  X  10-2 

4 

10.117  x  10-6 

2G 

5.2  X  10-Z 

5 

10.521  x  10" 6 

3G 

7.2  X  10-2 

6 

33.354  x  10-6 

Note  this 

is  a  steady  current. 
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FIELD  (1x10  G) 


18. 1985 

18.1585 

18.1185 

FIELD 

50.00 

49.99 

49.98 


FIG.  V- 
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SAMPLE  CELL  BOUNDARIES 
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SAMPLE  CELL  BOUNDARIES 


&  AXIS  COIL 


4  AXIAL  FIELD  PROFILES.  Calculated  field  values  along 
each  coil  pair  axis,  over  the  sample  cell  goundaries 


E.  EXPERIMENTAL  PROCEDURE 


The  basic  proced'— e  for  running  this  experiment  would  be  as  follows. 

1)  Turn  on  the  5b-axis  coils  and  condense  the  polarized  He3  - 

4 

He  mixture  in  the  sanple  cell.  At  this  point  the  magnetiza¬ 
tion  should  be  aligned  with  the  applied  field, 

2)  Turn  off  the  Sb-axis  coils  and  turn  on  the  ^G— axis  coils 
at  a  steady  value.  At  this  point  "iJ  should  be  precessing 
about  the  ib-axis. 

3)  Monitor  the  precessing  via  the  magnetometer  output. 

When  the  magnetometer  output  indicates  a  maximum  (or  minimum), 
turn  off  the  ^G-axis  coils  and  turn  on  the  3G-axis  coils. 

Now  should  be  precessing  about  the  ib-axis.  This  is  the 
normal  operating  condition  for  the  gyro. 

The  next  step  is  to  apply  the  rotating  field.  Figure  V-6  shows  a 
simple  method  of  generating  this  rotating  field  using  a  TR-20  desktop 
analog  computer.  To  establish  the  initial  conditions  described  in  the 
previous  chapter,  follow  this  procedure. 

1)  With  the  analog  computer  in  the  reset  mode,  monitor  the  magnet¬ 
ometer  output.  When  the  magnetometer  indicates  a  maximum, 
close  switch  No.  1. 

2)  Continue  to  monitor  the  magnetometer  output;  and,  at  the 
next  maximum,  switch  the  analog  computer  to  the  operate  mode. 

3)  To  change  the  direction  of  u  ,  toggle  switch  No.  2. 

The  peak  current  for  each  coil  is  determined  by  the  particular  combina¬ 
tion  of  initial  condition  (IC)  and  scale  factor  (SF)  pots. 
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FIG.  V-6  TR-20  CIRCUIT  TO  GENERATE  Bj. 

Potentiometers  cuz  set  the  coning  frequency;  Pot 
IC  sets  the  amplitude  and  is  used  in  conjunction 
with  pots  SF,  and  SF„  to  determine  peak  currents 
in  the  1G  and  2G  axis  field  coils. 
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Chapter  VI 


FIELD  GRADIENT  EFFECTS 

A..  FIELD  GRADIENT  SOURCES 

3 

Recall  that  the  basis  for  determining  rotations  with  a  He  nuclear 
gyro  is  the  stability  and  uniformity  of  the  applied  magnetic  field. 
Stability  is  enhanced  by  use  of  superconducting  field  coils  and  super¬ 
conducting  magnetic  shielding;  but  we  have  already  noted  that  our 
present  capabilities  do  not  permit  us  to  trap  a  perfectly  uniform  field. 
We  also  saw  in  the  previous  chapter  that  the  field  generated  by  a  helm- 
holtz  pair  is  not  perfectly  uniform.  This  latter  source  of  nonuniformity 
can  be  reduced  by  using  larger  helmholtz  coils  or  going  to  a  higher 
order  helmholtz  configuration  which  would  require  additional  coil  pairs 
to  correct  for  nonuniformity.  So  if  we  are  willing  to  pay  the  price 
on  the  field  coils,  we  are  limited  by  the  nonuniformity  of  the  trapped 
field.  Based  on  present  accomplishments  [Cabrera  1975]  we  would  expect 
to  be  able  to  obtain  fields  with  a  linear  gradient  term  of  less  than 
10_8G  cm_1  (10_10T  m_1). 

B.  REDUCTION  OF  TRANSVERSE  RELAXATION  TIME 

One  very  significant  effect  of  magnetic  field  gradients  i  t  t  - 
duce  the  transverse  relaxation  time,  T  .  Taber  [Taber  197' ^  n&s  co. .red 
the  theory  of  this  effect  in  detail,  so  only  the  final  result  and  a 
qualitative  explanation  are  included  here.  Consider  the  situation  shown 
in  Fig.  VI-1.  The  magnetic  field  has  a  linear  axial  gradient  such  that 

t  =  (BQ  +  az)z  (6.1) 

where  a  is  the  linear  gradient  coefficient.  Since  the  precession 
frequency  for  a  given  nucleus  depends  on  the  local  value  of  the  field, 

3 

the  two  He  nuclei  shown  in  Fig.  VI-1,  though  initially  aligned,  are  not 
processing  at  the  same  frequency.  If  these  nuclei  remained  in  their 
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X 


SAMPLE  CELL 


FIG.  VI-1  AXIAL  FIELD  PROFILE  ACROSS  SAMPLE  CELL. 
3  =  (BQ  +  az)z. 
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respective  positions,  then  the  angle  between  their  spin  axes  would  grow 

3 

linearly  with  time.  However,  the  He  nuclei  are  not  stationary,  but 
are  continuously  diffusing  through  the  sample  cell.  If,  on  the  average, 
a  given  nucleus  will  spend  as  much  time  in  one  location  as  any  other, 
then  the  gradient  is  averaged  out,  and  the  average  precession  rate  is 
TBp  rad/sec.  Though  each  will  have  the  same  average  precession  fre¬ 
quency,  the  two  nuclei  shown  in  Fig,  VI-1  will  have  different  magnetic 
field  time  histories  and  thus  will  still  become  misaligned.  The  faster 
they  diffuse  through  the  sample  cell,  though,  the  shorter  will  be  the 
time  required  to  average  out  any  differences,  and  thus  the  longer  will 
be  the  time  that  they  maintain  close  relative  alignment. 

3 

It  is  this  net  alignment  of  a  large  number  of  He  nuclei  that  pro¬ 
duces  the  signal  sensed  by  the  magnetometer.  So  as  the  alignment  becomes 
scrambled,  the  magnetization,  and  hence  the  signal  decays  to  zero.  The 
rate  of  this  decay  is  exponential  with  a  characteristic  time,  T 
given  by 

T~J  £  0.16Y2  a2  R4q  D-1  (6.2) 

3 

where  D  is  the  He  diffusion  coefficient.  So  to  maximize  T  we 

zu 

would  make  RQ  very  small  and  D  very  large.  However,  we  cannot  do 
this  with  impunity.  D  depends  on  density  and  on  temperature;  and  since 
the  field  coils,  pickup  coils  and  magnetic  shields  are  superconducting 
we  are  limited  in  temperature.  Further,  if  it  turns  out  to  be  advan¬ 
tageous  to  coat  the  sample  cell  walls  with  solid  hydrogen,  then  we  must 

maintain  the  temperature  below  8  K,  We  also  have  restrictions  on  the 
3 

He  density.  For  a  given  polarization  and  volume,  the  total  magnetic 
flux  produced  by  the  sample  is  proportional  to  its  density.  Thus  to 
increase  D  by  lowering  the  density  would  require  increasing  RQ  to 
maintain  the  same  signal  level.  So  one  must  trade  off  decay  rate  versus 
initial  magnitude.  Then,  as  we  will  see  in  the  next  two  sections,  there 
are  other  factors  which  will  influence  the  choices. 

Note  that  this  depolarization  process,  referred  to  as  gradient 
induced  relaxation,  affects  the  transverse  but  not  the  longitudinal 
relaxation  time.  Since  nuclei  aligned  with  the  field  (or  anti-aligned) 
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do  not  preeess,  their  alignment  is  not  affected  by  any  axial  field 
gradients.  This  is  one  factor  which  makes  the  transverse  relaxation 
time  generally  shorter  than  the  longitudinal  relaxation  time. 

£.  GRADIENT  INDUCED  g-SENSITIVE  DRIFT 

In  addition  to  shortening  the  relaxation  time,  magnetic  field 

gradients  have  a  more  subtle  influence  on  nuclear  gyro  performance. 

Should  the  external  environment  (e.g.,  thermal  or  gravitational)  induce 

3 

a  density  gradient  in  the  He  which  has  a  component  along  the  field 

gradient,  then  the  average  field  seen  by  the  nuclei  will  no  longer 

be  the  field,  B^,  at  the  geometric  center  of  the  sample  cell.  Thus 

there  will  be  a  shift  in  the  average  precession  rate,  which,  to  the 

extent  that  it  is  neither  measured  nor  corrected,  will  result  in  a 

drift.  These  effects  will  be  more  pronounced  in  gaseous  samples  but 

3  4 

will  also  be  present  in  liquid  He  -  He  solutions. 

So  let  us  consider  the  case  of  a  gaseous  sample  with  an  axial 
field  gradient,  undergoing  a  constant  acceleration  along  the  field  direc¬ 
tion.  In  particular,  consider  the  field  given  by  (6.1)  and  the  accel- 

^  3 

eration  to  be  in  the  z  direction.  Then  the  He  density  will  be  a 
function  of  z  and  can  be  determined  as  follows. 

Consider  the  one— dimensional  equation  of  hydrostatic  equilibrium 


dP 

dz 


Then  if  we  consider  the  He 


=  -Pg  . 

to  be  isothermal 


(6.3) 


and 


P  =  f&T 


££ 

P 


RT 


dz 


(6.4) 


(6.5) 


The  quantity  RT/g  is  commonly  referred  to  as  the  scale  height, 
Hq,  and  we  can  write 


d£ 

P 
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(6.6) 


If  we  integrate  over  the  sample  cell  from  z  =  -R  to  z,  we 
get  the  density  as  a  function  of  z  to  be 


P(z)  =  P(-RQ)e_(Z+Ro^Ho  .  (6.7) 

3 

Now  we  can  compute  the  average  magnetic  field  seen  by  the  He  nuclei 
as 


B 

av 


f pin  dv 
/  Pdv 


(6.8) 


where  the  integrations  are  carried  out  over  the  sample  cell  volume. 
Since  both  1?  and  p  are  functions  only  of  z  we  can  rewrite  (6.8) 
as 


B 

av 


°  P(-RQ)e  (Z+Ro)'/Ho  (bq  +  az)«(R2  -  z2)dz 

2% _ 

(R 0  -(z+R  )/H  „  2 

\  P(-R  )  e  n(R2  -  z  )dz 

— R  _  °  0 


(6.9 


Performing 


the  integrations  gives  the  results 


B  =  +  aR„ 

av  0  0 


— r  2  r  2 

e  (r  +3r+3)-e(r  -3r-t-3) 


r[e  r(l  +  r)  -  er(l  -  r)] 


where 


r 


R 

fT 


0 

0 


(6.10) 


(6.11) 


So  if  we  rewrite  (6.10)  as 


Bav  =  B0  +  58  (6-12) 

we  see  that  the  result  of  a  constant  acceleration  is  to  shift  the 
average  magnetic  field  by  5B.  This  will  cause  the  precession  fre¬ 
quency  to  shift  by  75B  rad/sec. 
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We  could  plot  families  of  curves  showing  the  drift  rate  vs  accel¬ 
eration  level  for  a  number  of  values  of  RQ,  a,  and  T,  but  it  is 
more  convenient  to  plot  a  single  curve  using  two  nondimensional  quan¬ 
tities.  These  are  a  normalized  drift  rate  YBB/YaR^,  which  represents 
a  fraction  of  the  maximum  drift  rate  due  to  acceleration  and  the  sample 
cell  radius  as  a  fraction  of  scale  height  (R^/H^).  This  one  curve 
(Fig.  VI-2)  contains  all  the  information  to  determine  the  g-sensitive 

drift  rate  for  any  set  of  conditions.  Figure  VI-2  shows  that  the 

-1  2 

drift  rate  increases  with  g,  T  ,  and  R 

As  an  example,  it  is  interesting  to  calculate  this  effect  for  the 
3 

He  nuclear  gyro  being  designed  for  the  electric  dipole  moment  experi- 

_  2 

ment.  Assuming  a  sample  cell  temperature  of  5K ,  R  of  1,9  x  10  m, 
lg,  and  a  linear  field  gradient  of  10  3G-cm  *  (10  10T-m  ^),  we 
get 

1,413  X  103  m 
1.345  x  10~5 

2.8  X  10-6  (from  Fig.  VI-2) 

—4 

S72  X  10  rad/sec  . 

_  q 

Thus  the  drift  rate  Y5B  is  1.084  X  10  rad/sec  (0.00022  deg/hr). 
Then,  since  this  effect  is  linear  in  acceleration,  we  get  for  this 
device  a  g-sensitive  drift  term  of  0.00022  deg/h  per  g.  And  this  is 
a  worst  case  where  both  the  field  gradient  and  the  density  gradient  are 
parallel.  In  general,  the  g-sensitive  drift  term  will  scale  as  the 
cosine  of  the  angle  between  the  two  gradients. 


Y5B_ 

7a  Fi 


YaR 
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_D.  GRADIENT  INDUCED  TEMPERATURE  SENSITIVE  DRIFT 

Now  let  us  evaluate  the  effect  of  a  temperature  gradient  induced 
density  gradient.  Again,  use  Eq.  (6.1)  for  the  magnetic  field  and 
describe  the  temperature  as 


T(z)  =  Tq  +  bz 


(6.13) 


where  T^  is  the  temperature  at  the  center  of  the  sample  cell  and 

/ 

b  is  a  linear  temperature  gradient  coefficient  (assumed  in  the  z-axis 
direction  to  give  a  worst  case  drift).  Now  setting  the  acceleration  to 
zero,  (6.3)  gives  a  constant  pressure.  Then  from  (6.7)  we  get 


where 


P(z) 


p  T 

0  0 


T  +  bz 
0 


po  ■  pb~1t01  • 


(6. 14) 


(6.15) 


So  again  with  both  and  p  as  functions  only  of  z 
as 


B 

a  v 


»<*o  -  *  ■*>  - 


R 


OPT 

-°  °  Jt (Rq  -  z  )  dz 

-Rq  T0  +  bZ  ° 


we  can  write  (6.8) 


(6.16) 


And  performing  the  integrations  gives 


B  +  aR 
0  0 


4  -1 

-Jt  -  2t  -in 


1 4  t 


-1 


l-t~ 


rr  d-t  ) 


2  +  In 


»)  - 


t) 


(6.17) 


where  t  is  the  temperature  at  the  center,  T  ,  divided  by  the  maxi¬ 
mum  temperature  difference  from  T  ,  bR^  (i.e.  t  =  T^/bR^). 
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Rewriting  (6.17)  as 


B  =  B„  +  gB’ 
av  0 


(6.18) 


it  is  again  convenient  to  plot  nondimensional  quantities  ygB'/yaR 

-1  0 
vs  t  (Fig,  VI-3)  rather  than  families  of  curves. 

As  an  example  let  us  again  use  the  gyro  design  of  the  electric 
dipole  moment  experiment  and  calculate  how  large  a  temperature  gradient 
can  be  allowed  while  keeping  the  temperature  gradient  induced  drift  less 

_9 

than  0.001  deg/hr  (4.85  X  10  rad/sec).  Using  the  same  values  as 
before,  we  get 


7aR0 

—4 

=  3,872  x  10  rad/sec 

78B’ 

TaR0 

=  1.252  X  10~5 

t"1 

bR 

=  — —  =  6.4  X  10  (from  Fig 

xo 

Ro 

To 

=  3.8  X  10-3  m-K-1  . 

Thus 

b 

=  1.68  X  10-1  K-nf1  . 

This  would  require  temperature  control  of  about  6 


VI-3) 


parts  in 


10 
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Chapter  VII 

CONCLUSIONS  AND  RECOMMENDATIONS 


A_.  CONCLUSIONS 

The  mathematical  models  developed  in  this  research  may  be  used  to 

3 

predict  the  behavior  of  a  cryogenic  He  nuclear  gyro  in  the  presence  of 
rotations,  vibrations,  and  magnetic  field  gradients  (i.3.,  field  gradient 
interactions  with  acceleration  ami  thermal  gradients). 

Application  of  these  models  to  the  electric  dipole  moment  gyro 
showed  that  thermal  gradient  control  requirements  would  not  be  stringent 
(~0.06%);  and  sensitivity  to  acceleration  is  very  small  (0.00022  deg/h-g). 

The  major  problems  uncovered  through  exercising  these  models  relate 
to  oscillatory  cross  axis  rotations.  These  types  of  motions  introduce 
drifts  through  various  rectification  mechanisms  and  can  lead  to  a  loss 
of  output  signal  by  causing  the  sample  magnetization  to  align  with  the 
applied  field.  The  latter  situation  can  be  precluded  through  the  use 
of  additional  field  coils  in  the  cross  axes  to  reduce  cross  axis  rota¬ 
tion  effects;  and,  if  necessary,  to  torque  the  sample  magnetization 
away  from  the  applied  field  direction. 

With  regard  to  the  rectification  mechanisms  we  find  conflicting 
requirements.  Drifts  due  to  high  frequency  vibrations,  that  reduce  the 
effective  field,  are  reduced  by  decreasing  the  applied  field.  Cross 
axis  rectification  effects,  however,  are  reduced  by  increasing  the 
applied  field  .  Also,  since  increasing  the  applied  field  increases 
the  precession  frequency  of  the  magnetization,  increasing  the  applied 
field  will  Increase  the  bandwidth  of  the  gyro  and  thus  increase  the 
effectiveness  of  the  cross  axis  field  coils.  One  solution  to  this 
apparent  dilemma  would  be  to  generate  as  large  an  applied  field  as 
practical  (field  gradients  generally  increase  with  field  magnitude) 
and  then  employ  a  mechanical  mount  which  would  attenuate  high  frequency 
cross  axis  oscillations.  The  mechanical  mount  should  be  notched  to 
severely  attenuate  vibrations  at  the  resonance  frequency. 
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The  experiment  discussed  in  Ch.  V  allows  verification  of  the  fre¬ 
quency  dependent  of  both  kinematic  and  cross  axis  rate  rectification 
mechanisms  by  varying  the  coning  frequency.  This  experiment  may  be 
conducted  with  a  relatively  simple  apparatus;  and,  by  using  rotating 
magnetic  fields,  no  mechanical  coning  motions  are  required.  The  smallest 
induced  drift  for  this  experiment  was  greater  than  2  deg/min.  This  is 
easily  observed  and  over  an  order  of  magnitude  greater  than  drifts  due 
to  other  sources. 

jB.  RECOMMENDATIONS 

There  are  a  number  of  areas  where  technological  improvements  will 
3 

make  the  cryogenic  He  nuclear  gyro  more  practical.  Among  these  would 

be  development  of  small  size  closed  cycle  cryogenic  refrigerators, 

3 

development  of  a  method  to  polarize  the  He  in  the  sample  cell  at 

cryogenic  temperature,  and  development  of  other  cryogenic  devices  such  as 

accelerometers,  gravity  gradiometers ,  and  computers  that  could  be 
3 

packaged  with  the  He  gyro  to  provide  an  all  cryogenic  inertial  measure¬ 
ment  unit.  Development  of  methods  to  improve  the  magnetic  field  qual¬ 
ities  and,  in  particular,  to  increase  the  magnitude  of  the  applied 
field  without  producing  large  gradients  will  lead  directly  to  improved 

3 

performance  of  a  He  gyro. 

In  addition  to  the  experiments  planned  for  both  the  Macor  and  fused 
quartz  gyros,  it  would  be  highly  desirable  to  perform  some  experiments 
to  observe  cross  axis  effects.  Effective  cross  axis  rate  inputs  could 
be  applied  as  fields  through  the  cross  axis  field  coils  or  by  using  the 
2(i-axis  coll  for  the  applied  field  and  rotating  the  dewar  about  the 
3”G  axis. 

Development  of  a  three  orthogonal  gyro  simulation  (preferably  with 
the  capability  of  substituting  the  actual  gyro  for  one  of  the  simulated 
gyros)  would  allow  evaluation  of  various  feedback  schemes  to  reduce 
cross  axis  effects,  as  discussed  in  Chapter  III. 

The  nuclear  gyro  shows  its  greatest  potential  for  a  practical  appli¬ 
cation  in  a  strapped-down  mode.  However,  before  this  potential  can  be 
realized  the  cross  axis  effects  must  be  removed. 


-76- 


REFERENCES 


Abraham,  A.  The  Principles  of  Nuclear  Magnetism,  Oxford  Press,  1973. 

Cabrera,  Bias,  "The  Use  of  Superconducting  Shields  for  Generating  Ultra- 
Low  Magnetic  Field  Regions  and  Several  Related  Experiments,"  Ph.D. 
Thesis,  Dept,  of  Physics,  Stanford  University,  Stanford,  CA.,  94305, 
1975. 

Cannon,  R.H.,  "Kinematic  Drift  of  Single-Axis  Gryoscopes,"  J.  of  Applied 
Mechanics .  Paper  No.  57-A-72. 

Cohen-Tannoud j i  Claude,  et  al.,  Quantum  Mechanics.  (Vols.  1  and  2), 

John  Wiley  &  Sons,  New  York,  NY,  1977. 

Colegrove,  F.D.,  et  al.,  Physics  Rev. ,  Vol .  132,  pp.  2561,  1963. 

Fairbank,  W.M. ,  and  Hamilton,  W.O. ,  Proc .  10th  Int.  Conf.  on  Low  Temp¬ 
erature  Phys . ,  Vol.  1IB,  p.  327,  1967. 

Giffard,  R.P.,  J.  of  Low  Temp.  Phys . ,  Vol.  6,  pp.  533,  1972. 

Kanegsberg,  E.,  et  al.,  "investigation  of  Bias  Stability  and  Cross  Axis 
Effects  and  Brassboard  Gyro  Conceptual  Design  for  a  Nuclear  Magnetic 
Resonance  Gyro,"  Technical  Rept.,  Air  Force  Avionics  Lab.,  AFAL/ 

RWA,  WPAFB ,  Ohio  45433,  No.  AFAL-TR- 78-1 155 ,  Oct.  1979. 

Macek,  W.M.,  et  al.,  "Ring  Laser  Rotation  Rate  Sensor,"  Proc,  of  the 
Symposium  on  Optical  Masers,  New  York,  1963. 

SHE  Manufacturer  Handbook,  SHE  Corporation,  San  Diego,  CA. ,  92121. 

Silver,  A.H.,  Zimmerman,  J.E.,  Phys .  Rev, ,  Vol.  157,  pp.  317,  May 
1967. 

Simpson,  J.H.,  et  al.,  "Nuclear  Magnetic  Resonance  Gyro  Design  and 
Test  Study,"  Technical  Rept.,  Air  Force  Avionics  Lab.,  AFAL/RWA, 
WPAFB,  Ohio,  45433,  No.  AFAL-TR- 79- 11 99,  Dec.  1979. 

Taber,  M.A.,  "I:  Spin-Lattice  Relaxation  of  Dilute  Solutions  of 

Polari led  He^  in  Liquid  He^  in  Low  Magnetic  Fields  at  4K.  II:  An 
Analysis  of  a  proposed  Cryogenic  He3  Nuclear  Gyroscope  and  its 
Application  to  a  Nuclear  Electric-Dipole  Moment  Experiment,"  Ph.D. 
Thesis,  Dept,  of  Physics,  Stanford  University,  Stanford,  CA, 

94305,  1978. 

Zimmerman,  J.E.,  et  al.,  J,  of  Applied  Phys.,  Vol.  41,  No.  4, 
pp.  1572-1580,  March  1970, 


-77- 


END 


DATE 

FILMED 


DTIC 


